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ABSTRACT:  This  report  summarizes  an  investigation  of  the  seismic  response  of  Corps  powerhouse 
structures.  Research  results  include  (1)  the  characterization  and  organization  of  the  range  in  geometry  of 
Coips  powerhouse  substructures  for  generator  bays  and  service/erection  bays;  (2)  development  of  repre¬ 
sentative,  idealized  substructure  geometries  for  these  types  of  structures;  (3)  construction  of  a  series  of 
finite  element  models  of  the  substructures  and  subjecting  these  substructures  to  earthquake  acceleration 
time-histories  representing  central/eastem  United  States  ground  motions  and  ground  motions  of  the  Pa¬ 
cific  Northwest  of  the  United  States;  (4)  development  of  a  database  of  first-mode  periods  of  the  typical 
powerhouse  substructures;  (5)  computation  of  top-of-substructure  acceleration  time-histories  and  conver¬ 
sion  of  these  time-histories  into  response  spectra;  (6)  development  of  amplification  factor  curves  (as  a 
function  of  frequency)  that  quantify  the  magnitude  of  the  amplification  of  ground  shaking  by  the  sub¬ 
structures;  (7)  development  of  a  methodology  for  performing  a  simplified  seismic  analysis  of  the  super¬ 
structure  (or  of  the  equipment  on  top  of  the  substructure)  that  takes  into  account,  in  an  approximate  sense, 
the  effect  of  the  substructure  on  the  response  of  the  superstructures  in  the  decoupled  seismic  structural 
analysis. 
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Introduction 


1.1.  Background  to  Powerhouse  Substructure 
Seismic  Amplification  Study 

Powerhouses  are  hydraulic  structures  used  to  generate  power.  A  powerhouse 
contains  numerous  types  of  electrical  and  mechanical  equipment  as  well  as  a 
single  bay  building  system  (referred  to  as  the  superstructure  in  this  report)  used 
to  protect  this  equipment.  The  performance  of  both  the  equipment  resting  on  top 
of  the  substructure  and  the  superstructure  during  and  after  seismic  shaking  is  of 
concern  to  the  U.S.  Army  Corps  of  Engineers. 

In  recent  decades  technology  in  the  fields  of  earthquake  engineering  and 
structural  dynamics  has  matured  into  viable  methodologies  that  can  be  applied  in 
a  practical  fashion  to  the  seismic  evaluation  of  powerhouse  superstructures  and 
their  equipment.  Regrettably,  most  of  the  79  powerhouses  for  which  the  Corps  is 
responsible  across  the  United  States  were  designed  and  constructed  prior  to  this 
point  in  time.  Concern  regarding  the  performance  of  powerhouse  superstructures 
and  its  contents  during  earthquakes  led  the  U.S.  Army  Corps  of  Engineers  under 
the  Corps  Earthquake  Hazards  Reduction  Program  (CEHRP)1  to  initiate  research 
leading  to  the  development  of  simplified  procedures  for  conducting  seismic 
response  evaluations  of  powerhouses  and  their  contents. 


1.1.1  Complete  dynamic  analysis 

A  standard  structural  dynamics  approach  to  the  seismic  response  evaluation 
of  a  structure  such  as  a  powerhouse  is  to  (a)  develop  a  detailed  finite  element 
model  of  the  powerhouse  substructure  (i.e.,  the  mass  concrete  portion  identified 
as  the  substructure  in  this  report),  (b)  include  in  the  finite  element  model  the 
superstructure,  and  (c)  subject  this  complete  model  to  earthquake  shaking.  Earth¬ 
quake  shaking  is  typically  represented  by  either  an  acceleration  time-history  or 
by  a  response  spectrum  in  the  dynamic  finite  element  analysis.  An  example  of  a 
seismic  evaluation  of  a  Corps  powerhouse  substructure  and  superstructure  in  a 
single  finite  element  model  (referred  to  as  a  complete  analysis  in  this  report)  is 
described  in  Ebeling  et  al.  (2002)  for  St.  Steven  Powerhouse.  A  modification 
could  easily  be  made  to  the  dynamic  finite  element  model  by  including  key 


1  Symbols  and  unusual  abbreviations  used  in  this  report  are  listed  and  defined  in  the  Notation, 
page  xviii. 
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equipment  on  top  of  the  substructure  in  the  finite  element  model  in  order  to 
predict  and  evaluate  its  response  to  earthquake  shaking. 

1.1.2  Decoupled  dynamic  analysis 

A  decoupled  seismic  analysis  may  be  undertaken  if  the  mass  and  stiffness  of 
the  superstructure  and/or  equipment  are  insignificant  compared  with  the  mass 
and  stiffness  of  the  substructure;  and  the  presence  or  absence  of  the  superstruc¬ 
ture  (and/or  equipment)  does  not  alter  the  dynamic  response  of  the  substructure. 
Experience  with  the  seismic  analysis  of  Corps  powerhouses  by  the  senior  author 
of  this  report  indicates  that  a  decoupled  structural  analysis  approach  is  a  viable 
method  of  analysis  for  the  Corps  powerhouses  provided  that  the  response  of  the 
substructure  is  properly  taken  into  account  in  the  simplified  analysis.  The  advan¬ 
tage  of  a  decoupled  analysis  is  that  modeling  efforts  are  concentrated  on  the 
superstructure  with  the  influence  of  substructure  response  accounted  for  in  an 
approximate  manner. 

1.2.  Research  Objectives 

The  purposes  of  this  research  are  to  investigate  and  characterize  the  seismic 
response  of  a  range  in  idealized  Corps  powerhouse  substructures  to  earthquake 
shaking  and  to  develop  a  simplified  method  for  accounting  for  the  effect  of  the 
substructure  in  a  decoupled  seismic  analysis  of  the  superstructure  and/or  equip¬ 
ment  resting  on  top  of  the  powerhouse  substructure.  The  following  were  specific 
objectives: 

•  Characterize  and  organize  the  range  in  geometry  of  Corps  powerhouse 
substructures  for  generator  bays  and  service/erection  bays. 

•  Develop  representative,  idealized  substructure  geometries  for  each  of  the 
two  general  categories. 

•  Construct  a  series  of  finite  element  models  of  the  substructures  and 
subject  these  substructures  to  earthquake  acceleration  time-histories 
representing  central/eastem  United  States  ground  motions  and  United 
States  Pacific  Northwest  ground  motions. 

•  Develop  a  database  of  first-mode  periods  of  the  typical  powerhouse 
substructures. 

•  Compute  top-of-substructure  acceleration  time-histories  and  convert 
these  time-histories  into  response  spectra. 

•  Develop  amplification  factor  curves  (as  a  function  of  frequency)  that 
quantify  the  magnitude  of  the  amplification  of  ground  shaking  by  the 
substructures. 

•  Develop  a  methodology  for  performing  a  simplified  seismic  analysis  of 
the  superstructure  (or  of  the  equipment  on  top  of  the  substructure)  that 
takes  into  account,  in  an  approximate  sense,  the  effect  of  the  substructure 
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on  the  response  of  the  superstructure  in  the  decoupled  seismic  structural 
analysis. 


1.3.  Complete  Structural  Dynamic  Analysis 

In  order  to  understand  what  constitutes  a  decoupled  seismic  analysis  of  the 
superstructure,  a  complete  dynamic  analysis  of  the  substructure  and  super¬ 
structure  in  a  single  finite  element  model  will  be  discussed  first.  A  schematic 
drawing  of  a  complete  dynamic  analysis  is  shown  in  Figure  1-1.  A  key  feature  of 
this  type  of  analysis  is  that  the  substructure  and  superstructure  (and/or  equip¬ 
ment)  are  idealized  in  a  single  model.  This  figure  shows  a  rock  response  spec¬ 
trum  used  as  input  at  the  base  of  the  substructure,  and  the  internal  forces  and 
displacements  are  calculated  within  the  superstructure  and/or  the  equipment.  The 
substructure  is  usually  modeled  using  continua  elements  and  the  superstructure 
(and/or  equipment)  as  multi-degree-of- freedom  (MDOF)  beam  elements  and 
lumped  masses.  Note  that  an  acceleration  time-history  may  also  be  used  to 
characterize  earthquake  shaking  in  a  seismic  analysis. 


1.4.  Decoupled  Structural  Dynamic  Analysis 

An  alternative  approach  for  the  dynamic  analysis  of  Corps  powerhouse 
superstructures  and/or  equipment  resting  on  top  of  the  substructure  is  a 
decoupled  analysis.  A  decoupled  dynamic  analysis  of  a  structural  system  consists 
of  decoupling  the  superstructure  and/or  equipment  from  the  substructure  and 
analyzing  each  separately.  First,  the  rock  ground  motion  (or  response  spectra)  is 
transferred  using  a  dynamic  analysis  from  the  base  of  the  model  to  the  top  of  the 
substructure  using  a  dynamic  response  model  of  the  substructure  (or  the  results 
from  a  dynamic  analysis  of  such  a  substructure  model).  In  this  research  a  time- 
history/modal  response  analysis  is  performed  on  each  substructure  model  and  a 
time-history  is  calculated  at  the  top  of  each  substructure.  The  calculated  time- 
history  is  converted  to  a  top-of-substructure  response  spectrum,  which  is  in  turn 
used  to  analyze  the  superstructure  and/or  equipment  model.  In  most  analyses  the 
acceleration  time-history  (or  response  spectrum)  at  the  top  of  the  substructure  is 
not  the  same  as  the  rock  acceleration  time-history  (or  response  spectrum)  due  to 
dynamic  response  effects  of  the  substructure.  This  top-of-substructure  response 
spectrum  is  then  used  as  input  to  the  simplified  finite  element  model  of  the 
superstructure  (or  equipment),  from  which  internal  structural  forces  and  displace¬ 
ments  are  computed. 

In  this  research  study,  a  series  of  time-history/modal  analyses  of  idealized 
powerhouse  substructures  was  conducted  to  compute  the  top-of-substructure 
acceleration  time-history  responses.  For  each  top-of-substructure  analysis,  a 
corresponding  acceleration  time-history  was  computed,  as  shown  in  Figure  l-2a. 
Using  the  Figure  l-2a  top-of-substructure  response  spectra  and  the  response 
spectra  for  rock,  an  amplification  factor  curve  versus  frequency  was  developed  at 
the  top  of  the  substructure,  as  shown  in  Figure  1  -2b.  Using  the  amplification 
factor  curve  times  a  user-defined  rock  response  spectrum,  an  “amplified”  (top-of- 
substructure)  response  spectrum  would  be  constructed,  as  shown  in  Figure  l-2c. 
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Figure  1-1 .  Complete  dynamic  analysis  of  substructure  and  superstructure 
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Figure  1-2.  Idealized  decoupled  dynamic  analysis  of  substructure  (Sheet  1  of  3) 


Chapter  1  Introduction 


cn 
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c.  Application  of  amplification  factor  curve  to  compute  top-of-substructure  response  spectra 
Figure  1-2. (Sheet  3  of  3) 


A  simplified  modal  analysis  of  the  superstructure  or  equipment  could  then  be 
conducted  using  this  amplified  response  spectrum.  By  using  this  amplified 
response  spectrum  as  the  earthquake  demand  in  the  finite  element  model  of  the 
superstructure,  the  simplified  analysis  would  approximately  account  for  the 
amplification  effects  of  the  substructure  during  earthquake  shaking  without  the 
need  for  directly  modeling  the  substructure  in  the  finite  element  model. 


1.5.  Evaluating  Nonlinear  Behavior  of 
Powerhouse  Superstructures  by  Coupled  or 
Decoupled  Response  Spectrum  Analysis 

Linear  elastic  response  spectrum  analyses  using  coupled  or  decoupled 
analytical  models  can  be  used  to  evaluate  the  performance  of  powerhouse 
superstructures  even  when  their  behavior  is  nonlinear.  The  performance 
evaluation  process  would  follow  Federal  Emergency  Management  Agency 
(FEMA)  356  (2000)  Linear  Dynamic  Procedure  (LDP)  analysis  procedures. 

In  LDP  analyses  a  response  spectrum  analysis  is  used  to  determine  the  force 
and  displacement  demands  on  powerhouse  superstructure  elements. 

Force-based  demand-to-capacity  evaluations  are  made  for  “brittle”  failure 
modes  (i.e.,  force-controlled  actions).  Force-controlled  actions  include  shear 
(diagonal  tension),  sliding  shear,  reinforcing  steel  anchorage,  and  fracturing  of 
tensile  reinforcement.  Performance  is  acceptable  if  the  ratio  of  force  demand  to 
force  capacity  (DCR)  is  less  than  one  (i.e.,  capacity  is  equal  to  or  greater  than 
demand).  It  is  recognized  in  the  LDP  analysis  process  that  capacity  may  decrease 
as  greater  displacement  ductility  demands  are  made  on  structural  elements. 

An  equal  displacement  response  or  equal  energy  response  is  used  to  estimate 
the  displacement  demands  on  powerhouse  superstructure  elements  that  behave 
inelastically.  Displacement  demand  (or  displacement  ductility  demand)  to  capac¬ 
ity  ratio  (DCR)  evaluations  are  made  for  “ductile”  failure  modes  (i.e.,  displace¬ 
ment-controlled  actions).  Flexure  is  the  displacement-controlled  action  of 
interest.  Performance  is  acceptable  if  the  DCRs  are  less  than  values  deemed 
acceptable  for  those  performance  levels  established  for  the  powerhouse  super¬ 
structure  evaluation. 

Force  and  displacement  demand  to  capacity  evaluations  were  conducted  as 
part  of  the  St.  Stephen  Hydropower  Plant  seismic  evaluation  (Ebeling  et  al. 
2002),  and  as  such  illustrate  the  application  of  an  LDP  type  analysis  procedure 
with  respect  to  performance  acceptability. 


1.6.  Time-History  Modal  Response  Analysis 

A  structural  linear  dynamic  analysis  was  performed  on  a  series  of  idealized 
powerhouse  substructures.  The  earthquake  loads  were  represented  by  accelera¬ 
tion  time-histories  at  rock  as  will  be  explained  in  Chapter  2  of  this  report.  A 
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time-history  modal  response  analysis  consists  of  solving  the  equation  of  motion 
(1-1)  using  numerical  methods: 

Mii(t)  +  Cii(t)  +  Ku(t)  =  -  MXg(t)  ( 1  - 1 ) 

In  Equation  1-1 ,  M,  C,  and  K  represent  the  mass,  damping,  and  stiffness 
matrices,  respectively,  of  an  MDOF  system  model.  For  these  dynamic  analyses,  a 
simplified  procedure  is  used  to  account  for  the  effect  of  the  upstream  and  down¬ 
stream  pool  of  water  in  the  dynamic  analyses;  Westergaard’s  hydrodynamic 
water  mass  is  added  to  the  structure  mass  in  the  mass  matrix  M.  In  Equation  1-1, 
ii(t)  ,  u(t)  ,  and  u(t)  represent  the  relative  acceleration,  velocity,  and  displace¬ 
ment  vectors,  respectively  (Headquarters,  U.S.  Army  Corps  of  Engineers 
(HQUSACE)  2003).  The  earthquake  ground  motion  acceleration  is  represented  in 
Equation  1  - 1  by  the  vector  Xg(t) . 

The  equation  of  motion  described  in  the  previous  paragraph  uses  step-by-step 
integration  solution  in  the  time  domain.  An  alternative  approach  found  in  tech¬ 
nical  literature  is  to  use  the  frequency  domain  procedure. 

A  time-history  of  acceleration  contains  acceleration  values  at  regular,  small, 
discrete  intervals  in  time.  In  many  instances,  the  maximum  acceleration  is  used  to 
quantify  the  level  of  dynamic  response.  A  graph  that  presents  the  maximum 
acceleration  for  different  structural  periods  (or  frequencies)  is  referred  to  as  a 
response  spectrum  (refer  to  Ebeling  1992).  Figure  1-3  summarizes  the  procedure 
for  the  development  of  a  response  spectrum  curve.  A  response  spectrum  curve  is 
developed  assuming  the  damping  value  of  a  structure.  If  total  time-history  accel¬ 
eration  is  used,  the  acceleration  in  the  response  spectrum  curve  is  referred  to  as 
total  acceleration  (SA)  (sometimes  referred  as  absolute  acceleration  (Ebeling 
1992)  and  spectral  acceleration  (Chopra  1995)).  When  relative  time-history 
acceleration  is  used,  the  response  spectrum  acceleration  is  called  pseudo-spectral 
acceleration  (PSA)  or  SA.  Table  1-1  lists  the  relationship  between  pseudo-spectral 
values  and  total  values.  The  difference  between  PSA  and  SA  is  that  the  SA 
values  take  into  account  the  damping  characteristic  of  the  structure,  while  PSA 
values  do  not.  Chopra  ( 1 995)  has  shown  that  the  difference  between  these  two 
values  is  small  for  short-period  systems  and  significant  only  in  long-period 
structures  with  high  damping  values. 

The  dynamic  analysis  computations  summarized  in  this  report  were  made 
using  the  structural  analysis  software  SAP2000  (Computers  and  Structures,  Inc. 
(CSI),  1997).  SAP2000  uses  the  mode  supeiposition  method  to  solve  the 
dynamic  equilibrium  equation  of  motion.  For  linear  transient  analysis,  the  soft¬ 
ware  assumes  a  linear  variation  of  the  time  function  and  the  response  is  computed 
by  closed  form  integration  of  the  modal  equation. 
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Figure  1-3.  Development  of  a  response  spectrum  curve  (from  Ebeling  1992) 


1.7.  Monoliths  Analyzed  (Generator  and 
Erection/Service  Bay  Monoliths) 

A  powerhouse  consists  of  an  erection/service  bay  monolith  and  a  series  of 
generator  bay  monoliths.  Most  of  the  storage  rooms  are  located  in  the  erection/ 
service  bay  monolith.  The  erection/service  bay  monolith  comprises  mass  con¬ 
crete  and  internal  frames  for  the  floors  of  rooms  and  partitions.  The  turbine 
generators  and  other  electrical  and  mechanical  equipment  are  located  in  the 
generator  bay.  Construction  drawings  for  54  Corps  powerhouses  were  reviewed. 
Appendix  A  lists  these  powerhouses  and  summarizes  their  key  geometry 
characteristics. 
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Table  1-1 

Definitions  of  Earthquake  Response  Spectrum  Terms  (from  Ebeling, 
1992) 


Svmbolc 

Definition 

Dcscrlntion 

S„  -  SD 

1  x(t)  1 ... 

Relative  displacenent  response  spectrum 
or  spectral  displacement 

SV 

1  *(t)  1 

Relative  velocity  response  spectrum 

SA 

^  ^  BAT 

Absolute  acceleration  response  spectrum 

Sv  -  PSV 

..  c  2s  e 

w  Sq  •  j  Sp 

Spectral  pseudo- velocity 

SA  -  PSA 

s„  -  s0 

Spectral  pseudo-acceleration 

w  is  th«  circular  frequency  of  vibration  of  the  undaaped  SDO?  systea  in  units 
of  radians  per  second. 

T  is  the  natural.  (or  fundamental)  period  of  vibration  of  the  undamped  SDOF 
system  in  units  of  seconds 

*-WT 


A  generator  bay  monolith  was  categorized  into  one  of  four  groups,  depend¬ 
ing  on  its  height.  For  each  group,  an  average  height,  an  upstream/downstream 
dimension  (transverse  to  the  axis  of  the  powerhouse),  a  plan  width  (along  the  axis 
of  the  powerhouse),  and  pool  and  tailrace  water  elevation  were  determined. 

Figure  1  -4  shows  a  typical  cross  section  of  a  generator  bay,  this  one  for  the 
Cheatham  Powerhouse.  The  plan  view  of  the  Cheatham  Powerhouse  generator 
bay  is  shown  in  Figure  1-5. 

Figure  1-6  shows  an  example  of  a  cross  section  of  an  erection/service  bay 
monolith,  this  one  for  the  Cordell  Hull  powerhouse.  Four  groups  of  erection/ 
service  bay  monoliths  were  created,  and  each  of  the  54  Corps  powerhouses  was 
assigned  to  one  of  the  groups,  depending  on  its  height. 


1.8.  Idealized  Two-Dimensional  Finite  Element 
Models 

In  this  study,  the  amplified  response  spectrum  at  the  top  of  the  substructure 
will  be  determined  in  order  to  construct  amplification  factor  curves  representative 
for  each  model.  In  order  to  perform  the  dynamic  analysis,  the  substructures  were 
idealized  into  two-dimensional  (2-D)  representative  finite  element  models. 
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Figure  1-4.  Cheatham  Powerhouse  generator  bay  cross  section 


A  2-D  composite  finite  element  model  was  selected  for  the  idealization  of  the 
generator  bay  substructure.  It  consists  of  five  structural  members  that  form  a 
three-dimensional  (3-D)  “box.”  In  Figure  1-7,  the  concept  of  the  composite 
model  is  presented.  The  walls  that  undergo  upstream/downstream  direction 
shearing  deformations  during  earthquake  shaking  are  referred  to  as  shear  walls 
(Walls  1  and  2).  The  base  slab  is  referred  to  as  3.  Walls  4  and  5  are  flexural  walls 
defining  the  upstream  face  and  the  downstream  face  of  the  generator  bay  mono¬ 
lith,  respectively.  The  analyses  will  be  a  simplified,  two-dimensional  analysis 
using  plane  strain  elements  subjected  to  an  upstream/downstream  ground  motion. 
Due  to  the  nature  of  the  simplified  2-D  model  of  a  single  generator  bay  monolith, 
shear  Wall  2  is  superimposed  on  shear  Wall  1  in  the  “composite”  model.  The 
nodal  degrees  of  freedom  for  the  “composite”  element  structure,  consisting  of  the 
five  structural  features,  are  idealized  in  this  figure. 

Two  simplified  2-D  finite  element  models  were  created  to  represent  the 
geometry  and  properties  of  the  erection/service  bay  substructure,  a  block  model 
and  a  block-frame-shear  wall  model.  Block-frame-shear  wall  models  have 
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Figure  1-5.  Cheatham  Powerhouse  plan  view  at  generator  area 
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Figure  1-6.  Cordell  Hull  Powerhouse  erection/service  bay  cross  section 


Figure  1-7.  The  idealized  two-dimensional  composite  finite  element  model  for  the 
generator  bay  substructure  used  in  this  study 


frame  elements  within  the  substructure  to  represent  the  partitions  inside  the 
powerhouse,  providing  flexibility  to  the  models  compared  to  the  block  models. 
For  the  three  models,  nine-node  plane  stress  elements  were  used  in  the  finite 
element  models.  The  models  were  created  using  SAP2000  software. 
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1.9.  Report  Content 


This  report  summarizes  the  results  obtained  from  a  series  of  time-history 
analyses  of  generator  bay  and  erection/service  bay  substructures.  The  develop¬ 
ment  of  the  earthquake  ground  motions  and  their  characteristics  are  presented  in 
Chapter  2.  Chapter  3  discusses  the  development  of  the  amplification  factor 
curves  for  the  generator  bay  substructure.  The  amplification  factor  curves  for  the 
block  and  block- frame-shear  wall  model  of  the  erection/service  bay  are  presented 
in  Chapter  4.  The  recommended  amplification  factor  curves  are  summarized  in 
Chapter  5  followed  by  two  examples;  (1)  a  comparison  example  between  a  com¬ 
plete  and  a  decoupled  analysis  and  (2)  a  demonstration  example  of  an  amplified 
rock  response  spectrum  developed  using  the  NEHRP  (2000)  procedure,  to  be 
used  to  analyze  a  simplified  frame  superstructure  model. 

Appendix  A  discusses  geometry  characteristics  for  54  of  the  Corps  power¬ 
houses.  Appendix  B  discusses  the  methodology  used  to  idealize  the  upstream 
substructure  wall  (4)  in  the  composite  model  of  the  generator  bay.  Appendix  C 
discusses  the  results  from  a  parametric  study  of  planes  stress,  plane  strain  and 
shell  elements.  Appendix  D  describes  the  computation  of  the  hydrodynamic 
added  mass.  Appendix  E  discusses  PSA  and  SA  curves.  Appendix  F  summarizes 
amplification  curves  for  2  percent  damping. 
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2  Development  of 

Earthquake  Time-Histories 


2.1.  Introduction 

Two  earthquake  load  cases  were  developed  for  this  study  using  design 
earthquake  ground  motions  for  a  powerhouse  site  in  the  Pacific  Northwest  of  the 
United  States  (PNW)  and  a  powerhouse  site  in  the  central-eastern  United  States 
(CEUS).  Each  case  was  selected  from  an  initial  suite  of  sites  considered  for  the 
respective  cases.  The  Pacific  Northwest  site  was  selected  from  the  suite 
Bonneville  (symbol  BnV),  John  Day  (JDay),  the  Dalles  (Dal),  Chief  Joseph  (CJ), 
and  McNary  (McN)  Dams.  A  central-eastern  U.S.  site  was  selected  from  the 
eastern  sites  St.  Steven  (StS),  Carters  Main  Dam  (CMD),  and  Richard  B.  Russell 
(RBR)  and  central  U.S.  sites  Barkley  (Brk),  Greers  Ferry  (GrF),  and  Bull  Shoals 
(BuS).  Site  parameters  (location  and  peak  ground  acceleration  (PGA)  for  a  range 
of  exposure  periods)  for  the  two  suites  of  sites  are  listed  in  Table  2-1  and  mapped 
in  Figures  2-1  and  2-2  for  the  PNW  and  CEUS  sites,  respectively.  These  figures 
also  display  the  historical  seismicity  and  any  nearby  active  faults  for  the  region  to 
provide  a  general  view  of  the  seismic  hazards  that  affect  these  sites.  The  final 
earthquake  loading  was  provided  as  acceleration  time-histories.  These  corre¬ 
sponding  earthquake  time-histories  were  based  on  probabilistic,  mean  level, 
equal  hazard,  5  percent  damped  acceleration  response  spectra  corresponding  to 
the  seismic  hazard  of  the  selected  site  for  an  exposure  period  of  1,000  years.  The 
acceleration  time-histories  were  developed  synthetically  using  these  initial 
response  spectra  as  a  target.  Each  resulting  time-history  was  further  modified  to 
provide  a  close  spectral  match  with  the  basis  target  spectra  over  the  period  range 
of  interest  for  the  structural  response  analyses. 


2.2.  Spectra  Selection 

This  section  describes  the  investigation  and  selection  of  the  target  spectra 
that  were  used  to  develop  two  “matching”  acceleration  time-histories  (i.e., 
spectrum  compatible),  one  each  for  the  PNW  and  CEUS  sites.  The  approach  was 
to  develop  site-specific  equal  hazard  spectra  that  are  more  suited  as  a  basis  to 
develop  acceleration  time-histories.  These  spectra  were  based  on  the  seismic 
hazard  curves  from  the  probabilistic  seismic  hazard  analysis  (PSHA)  conducted 
by  the  U.S.  Geological  Survey  (USGS)  for  development  of  the  2002 
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Table  2-1 

Powerhouse  Sites  and  Probabilistic  PGA  for  Return  Periods  500, 
1,000,  and  2,475  Years 

Powerhouse  Site 

Latitude 

Longitude 

PGA,  g 

Return  Period,  years 

500 

1,000 

2,475 

PNW 

Bonneville 

45.642 

-121.933 

0.1335 

0.1808 

0.2609 

John  Day 

45.717 

-120.685 

0.0935 

0.1283 

0.1939 

The  Dalles 

45.615 

-121.138 

0.0987 

0.1337 

0.1989 

Chief  Joseph 

47.996 

-119.627 

0.0862 

0.1222 

0.1918 

McNary 

45.930 

-119.295 

0.0817 

0.1162 

0.1851 

CEUS 

Eastern 

St.  Steven 

33.425 

-79.933 

0.2158 

0.5748 

1.2018 

Carters  Main  Dam 

34.613 

-84.685 

0.0705 

0.1177 

0.2273 

Richard  B.  Russell 

34.025 

-82.592 

0.0615 

0.0991 

0.1787 

Central 

Barkley 

33.022 

-88.220 

0.1965 

0.3230 

0.5489 

Greers  Ferry 

35.525 

-92.000 

0.1077 

0.1897 

0.3202 

Bull  Shoals 

36.363 

-92.573 

0.0631 

0.1063 

0.1756 

U.S.  National  Seismic  Hazard  Maps  (Frankel  et  al.  1996,  2000,  2002).  An 
alternative  would  have  been  to  use  generalized  design  spectra  such  as  provided  in 
FEMA  302  (FEMA  1997)  and  ASCE  7-2  (American  Society  of  Civil  Engineers 
2002)  guidance  for  seismic  analyses. 


2.2.1  Pacific  Northwest  site 

To  facilitate  selection  of  a  single  response  spectrum  to  represent  this  case, 

5  percent  damped,  equal  hazard,  pseudo-spectral  acceleration  response  spectra 
(EHPSA)  were  developed  for  a  suite  of  sites  (shown  in  Figure  2-1)  for  return 
periods  of  500,  1,000,  and  2,475  years.  It  was  decided  to  use  1,000-year  design 
ground  motions  as  the  more  appropriate  exposure  level  for  performance  objec¬ 
tives  considered  important  for  this  study.  Inspection  of  these  EHPSA  plotted  in 
Figure  2-3  shows  a  common  spectral  shape  anchored  to  range  of  PGA.  For  the 
1,000-year  return  period  the  PGA  ranges  from  0.1 16  to  0.180  g.  These  spectra 
peak  in  the  0.2-sec-period  range  with  a  broader  dispersion  of  0.27  to  0.44  g.  The 
spectra  for  the  Bonneville  site  is  an  upper  bound  for  the  1,000-year  spectra  and 
additionally  a  representative  lower  bound  to  the  2,475-year  spectra.  Further 
analysis  of  this  suite  of  spectra  is  presented  in  Figure  2-4  where  these  spectra  are 
replotted  as  a  ratio  of  PSA  and  PGA.  This  plot  shows  that  the  spectral  shape  is 
common  and  its  relationship  to  the  anchor  PGA  is  very  similar.  Again,  the 
selected  Bonneville  spectra  are  shown  to  be  an  upper  bound  in  the  shorter  period 
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Figure  2-1 .  PNW  powerhouse  sites  and  active  faults  and  seismicity  M  >  4.0 

(higher  frequency)  band,  although  not  anomalously.  For  this  study  the  higher 
frequency  components  of  this  loading  are  the  more  important.  This  plot  also 
shows  reasonableness  in  using  PGA  to  scale  ground  motions  for  different  ranges 
of  loads  while  maintaining  an  analogous  relative  frequency  content  for  this 
region. 

Also  presented  in  these  two  figures  for  comparison  are  the  selected  spectra 
for  the  CEUS  site.  This  comparison  shows  that  the  spectra  for  the  CEUS  site  are 
richer  in  high  frequency  as  evidenced  by  the  higher  PGA  and  shift  of  the  pre¬ 
dominant  frequencies  to  a  higher  band.  The  converse  can  be  seen  in  the  figures 
showing  the  spectra  for  the  CEUS  site. 

One  important  aspect  of  a  PSFIA  is  that  the  resulting  hazard  is  a  sum  of  all 
significant  sources  that  affect  a  given  location.  This  is  shown  in  Figure  2-5  as  a 
graph  of  deaggregated  seismic  hazard  for  the  selected  Bonneville  site.  The 
presented  component  of  hazard  for  deaggregation  was  the  0.2-sec,  974-year 
(nominal  1,000-year)  data.  It  can  be  seen  that  this  component  of  hazard  for  this 
site  is  composed  mainly  of  large-magnitude  M  =  8+  earthquakes  at  ranges 
>160  km  and  a  cluster  of  small-magnitude  M  <  7  earthquakes  within  80  km  of 
the  site.  These  two  distinct  clusters  of  sources  correspond  to  large  earthquakes 
from  the  Cascadia  subduction  zone  and  closer  shallow  crustal  faults,  respectively. 
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Figure  2-2.  CEUS  powerhouse  sites  and  seismicity  M  >  4.0 

The  calculated  mean  hazard  for  this  spectral  period  and  exposure  period  from  all 
these  sources  is  an  M  =  6.77  earthquake  at  a  distance  of  61  km.  This  would  not 
necessarily  be  the  same  for  all  spectral  periods  or  exposure  periods,  which 
explains  showing  the  component  most  appropriate  for  this  study:  1,000-year, 
spectral  peak  in  medium- high  frequency  range.  A  final  related  point  follows: 
selecting  a  single  recorded  acceleration  time-history  (based  on  a  single  magni¬ 
tude  and  distance)  to  represent  spectra  derived  from  a  summation  of  sources  at 
different  magnitudes  and  ranges  may  be  problematic.  Therefore,  using  an 
accelerogram  developed  synthetically  to  realistically  match  an  EHPSA  is  more 
reasonable,  especially  if  phasing  of  the  arriving  waves  is  not  an  important  aspect 
in  the  analysis,  which  is  the  case  for  this  study. 


2.2.2  Central-eastern  U.S.  site 

As  stated  previously,  to  facilitate  selection  of  a  single  response  spectrum  to 
represent  this  case,  5  percent  damped  EHPSA  were  developed  for  the  suite  of 
sites  mapped  in  Figure  2-2  for  return  periods  of  500,  1,000,  and  2,475  years.  The 
1,000-year  ground  motion  was  chosen  as  the  appropriate  exposure  level  for  the 
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Figure  2-3.  Northwest  powerhouse  site  equal  hazard  spectra 
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Figure  2-5.  Bonneville  site  probabilistic  hazard  deaggregation  for  5-Hz  spectral 
component 

performance  objective  considered  important  for  this  study.  These  EHPSA  plotted 
in  Figure  2-6  again  show  a  common  spectral  shape  anchored  to  range  of  PGA. 
For  the  1,000-year  return  period  the  PGA  ranged  from  0.099  to  0.574  g.  These 
spectra  peak  in  the  0.1 -sec-period  range  with  a  broader  dispersion  of  ~0.2  to 
1.25  g.  This  is  a  much  wider  spread  with  a  peak  at  a  higher  frequency  than  the 
selected  PNW  site.  The  Greers  Ferry  site  spectrum  was  selected  to  represent  the 
CEF1S  region.  Figure  2-7  shows  it  is  not  an  upper  bound  for  the  1,000-year 
spectra,  surpassed  by  St.  Steven  and  Barkley  sites.  However,  St.  Steven  and 
Barkley  were  considered  too  severe  for  a  general  CEUS  case.  The  PSA/PGA 
spectra  plots  in  Figure  2-7  show  that  the  Greers  Ferry  spectrum  has  a  relatively 
high  bounding  ratio  in  the  suite,  which  supports  it  as  a  good  choice,  especially 
for  scaling  load  cases  using  PGA. 

The  seismic  hazard  deaggregation  is  shown  in  Figure  2-8  for  the  selected 
Greers  Ferry  site.  The  presented  component  of  hazard  for  deaggregation  was  the 
0.2-sec,  974-year  (nominal  1,000-year)  data.  It  can  be  seen  that  this  component 
of  hazard  for  this  site  is  composed  mainly  of  large-magnitude  M  =  7.5+  earth¬ 
quakes  at  ranges  >150  km  and  a  broad  cluster  of  small-magnitude  M  <  7  earth¬ 
quakes  within  50  km  of  the  site.  These  two  clusters  of  sources  correspond  to 
large  earthquakes  from  the  New  Madrid  zone  and  closer  shallow  crustal  faults 
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Figure  2-6.  CEUS  powerhouse  site  equal  hazard  spectra 


Figure  2-7.  CEUS  powerhouse  site  spectral  ratios 
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Figure  2-8.  Greers  Ferry  site  probabilistic  hazard  deaggregation  for  5-Hz  spectral 
component 


associated  with  diffuse  seismicity,  respectively.  The  calculated  mean  hazard  for 
this  spectral  period  and  exposure  period  from  all  these  sources  is  an  M  =  7.26 
earthquake  at  a  distance  of  1 07  km. 


2.3.  Time-History  Development 

The  acceleration  time-histories  were  developed  in  a  two-step  process  once  a 
basis  or  target  site  and  EHPSA  were  selected:  (a)  develop  a  synthetic  seismogram 
constrained  to  the  site  EHPSA  and  (b)  further  modify  the  synthetic  accelerogram 
by  close  spectral  matching  with  the  site  EHPSA. 


2.3.1  Synthetic  development 

Considering  the  objectives,  scope,  and  type  of  dynamic  analysis  of  this  over¬ 
all  study,  it  was  decided  to  develop  the  time-histories  synthetically  instead  of 
using  recorded  accelerograms.  The  scope  was  to  employ  one  earthquake  time- 
history  load  for  the  site.  The  criteria  for  developing  the  load  case  was  based  on 
ground  motion  levels  with  a  5  percent  probability  of  exceedance  (PE)  for  an 
exposure  period  of  50  years  or  return  period  of  974  years,  which  was  set 
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nominally  at  1,000  years.  These  probabilistically  derived  ground  motion  levels 
were  specified  as  PGA  and  spectra  accelerations  as  uniform  or  equal  hazard 
response  spectra  for  a  mean  level  with  5  percent  damping.  Since  these  derived 
target  spectra  are  smooth  and  derived  from  summing  from  multiple  sources,  it 
was  considered  that  starting  with  a  deterministic  single  earthquake  source,  single 
recording,  would  require  considerable  spectral  modification.  Since  the  original 
character  of  a  recorded  accelerogram  would  be  basically  overwritten  to  get  the 
required  good  match,  this  negated  one  advantage  for  using  a  natural  record. 
Finally,  as  the  load  case  was  not  as  severe  as  a  Maximum  Credible  Earthquake 
(MCE)  level  event,  an  elastic  response  was  expected  with  the  results  not  as 
sensitive  to  the  phasing  of  the  load  as  a  nonlinear,  possibly  large  deformation 
response  from  an  extreme  event.  The  synthetic  seismograms  were  developed 
from  the  strong  ground  motion  model  SMSIM  TD  (Boore  2000).  This  is  a 
simple  model  based  on  a  stochastic  point  source,  which  is  reasonable  for  the 
constraints  and  objectives  of  this  study.  The  synthetic  seismogram  was  devel¬ 
oped  using  an  attenuation  relationship  and  source  spectra  appropriate  for  the 
region  of  the  site  scaled  (constrained)  by  the  target  spectra  associated  with  the 
spectral  component  and  probability  of  exceedance  of  the  EHPSA  of  the  selected 
site  (0.2-sec  PSA,  5  percent  PE  in  50  years). 


2.3.2  Spectral  Matching 

Spectral  matching  was  accomplished  using  the  RSPMATCH  computer 
program  (Abrahamson  1998)  that  is  based  on  the  time  domain  method  of  Tseng 
and  Lilhanand  (Lilhanand  and  Tseng  1988).  The  resulting  spectrally  matched 
synthetic  ground  motion  time-histories  are  shown  in  Figures  2-9  and  2-10  for  the 
PNW  and  CEUS  sites,  respectively.  Note  that  these  time-histories  are  generated 
at  0.005  sec,  which  can  provide  signal  fidelity  up  to  a  frequency  content  of  100 
Flz.  Flowever,  these  simulations  have  limited  frequency  content  (insignificant) 
above  30  Flz,  which  can  be  seen  by  inspection  of  the  frequency  versus  amplitude 
plots  shown  in  Figures  2-9  and  2-10.  Also  as  seen  in  these  figures,  the  time- 
histories  are  baseline  corrected  as  evidenced  by  the  zero  velocity  at  the  end  of  the 
time  series.  This  is  an  important  characteristic  that  avoids  errors  in  the  numerical 
analyses  due  to  inadvertent  long-period  drift  artifacts  in  the  accelerogram.  A 
more  detailed  plot  of  the  final  spectral  matches  showing  the  initial  synthetics 
from  SMSIM,  the  spectral  targets,  and  the  final  spectral  matched  spectra  for 
2  percent  and  5  percent  damping  are  shown  in  Figures  2-11  and  2-12  for  the 
PNW  and  CEUS  sites,  respectively.  Each  figure  contains  a  plot  of  four  5  percent 
damped  response  spectra: 

•  Seed  TH  5%  -  spectra  of  original  (or  seed)  synthetic  time-history  that 
was  to  be  modified  to  match  the  spectral  target. 

•  BnV  Target  5  percent  -  the  Bonneville  (BnV,  Figure  2. 1 1)  or  Greers 
Ferry  (GRF,  Figure  2-12)  site  target  spectra.  The  plotted  points  on  these 
target  spectra  were  the  specific  frequencies  that  controlled  the  spectral 
modification.  This  range  for  spectral  matching  was  from  0.03  to  2.0  sec 
(a  minimum  to  maximum  frequency  of  0.05  to  33  Hz). 
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Figure  2-9.  Acceleration  time-history  for  PNW  powerhouse  -  Bonneville  1 ,000-year  equal  hazard  spectra 
synthetic  and  spectrally  matched 

•  BnV  RSPMatch  2%  and  BnV  RSPMatch  5%  or  GRF  RSPMatch  2%  and 
GRF  RSPMatch  5%-the  resulting  2  percent  and  5  percent  damped 
spectra  from  the  accelerograms  spectrally  modified  to  approximately 
match  the  target  spectra. 
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Figure  2-10.  Acceleration  time-history  for  CEUS  powerhouse  -  Greers  Ferry  1 ,000-year  equal  hazard 
spectra  synthetic  and  spectrally  matched 
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Figure  2-1 1 .  Comparison  of  original,  target,  and  spectral  matched  spectra  for  the  PNW  powerhouse  site 
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Figure  2-12.  Comparison  of  original,  target,  and  spectral  matched  spectra  for  CEUS  powerhouse  site 
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3  Time-History  Finite 
Element  Results  for  the 
Generator  Bay  Monoliths 


3.1.  Introduction 

This  chapter  discusses  the  results  from  analyses  of  a  series  of  time-history/ 
modal  responses  of  the  four  categories  of  generator  bay  monoliths  to  earthquake 
loading.  The  earthquake  loading  is  represented  by  acceleration  time-histories  in 
the  analyses  discussed  in  this  chapter.  Two  acceleration  time-histories  were  used 
in  the  analysis,  one  representing  a  central-eastern  United  States  earthquake  event 
and  a  second  representing  a  Pacific  Northwest  event.  Each  acceleration  time- 
history  was  specified  at  the  base  of  the  finite  element  model  of  each  idealized 
powerhouse  substructure  (i.e.,  in-rock).  Two-dimensional,  plane  stress  finite 
elements  were  used  in  the  SAP2000  (CSI  1997)  finite  element  models  of  each  of 
the  four  idealized  generator  bay  models  studied  during  the  course  of  this 
research.  From  these  computed  time-history  results,  response  spectra  at  top-of- 
substructure  were  then  computed  and  converted  into  response  spectra.  These 
results,  in  turn,  were  converted  into  a  plot  of  response  spectra  amplification 
factors  (AF)  versus  frequency  for  use  in  decoupled  response  spectra  modal 
analyses.  The  AF  curves  presented  in  this  chapter  are  plotted  up  to  50  Flz.  Flow- 
ever,  the  ground  motion  used  in  this  study  has  frequencies  up  to  33  Flz.  The  use 
of  AF  curves  above  33  Flz  is  not  recommended. 


3.2.  Geometry  of  Four  Height  Monoliths 

The  U.S.  Army  Corps  of  Engineers  is  responsible  for  approximately 
79  powerhouses  across  the  United  States.  Drawings  for  54  of  the  generator  bays 
at  these  powerhouses  were  studied  and  assigned  to  one  of  four  groups  according 
to  height  of  substructure.  For  each  group,  average  upstream/downstream  widths 
(transverse  to  the  axis  of  the  powerhouse)  and  average  monoliths  widths  (parallel 
to  the  axis  of  the  powerhouse)  were  determined.  Additionally,  for  each  group,  a 
representative  normal  pool  elevation  (upstream  of  the  powerhouse)  and  a  repre¬ 
sentative  tailrace  elevation  (downstream  of  the  powerhouse)  were  determined. 
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The  pool  water  elevation  was  taken  as  a  constant  20  ft1,2  above  the  top-of- 
substructure  for  all  groups,  and  the  tailrace  water  elevation  varied  as  a  function  of 
group  height.  Table  3-1  lists  the  four  groups  with  the  values  assigned  and  their 
respective  pool  and  tailrace  water  elevations. 


Table  3-1 

Dimensions  of  Generator  Bay  Section 

Group 

Height  H 
ft 

Upstream/ 
Downstream 
Width  W,  ft 

Plan 

Width  B 
ft 

Pool 

Water  El 
ft 

Tailrace 

El,  ft 

Range  in  Height 
for  Group,  ft 

1 

125 

135 

85 

145 

89 

115-135 

2 

100 

130 

75 

120 

74 

90-115 

3 

70 

110 

60 

90 

60 

60-75 

4 

40 

65 

50 

60 

36 

25-60 

Appendix  A  contains  a  summary  table  of  the  dimensions,  pool  and  tailrace 
elevations,  and  general  characteristics  of  each  powerhouse  reviewed  in  this  study. 


3.3.  Material  Properties 

The  generator  bay  monoliths  are  massive  concrete  structures.  The  material 
properties  assumed  in  this  study  are  listed  in  Table  3-2  for  the  massive  concrete 
of  the  substructure.  Most  of  this  information  was  obtained  from  Ebeling  et  al. 
(2002). 


Table  3-2 

Material  Properties  Used  in  the  Generator  Bay  Mass  Concrete 

Substructure 

Mass  Concrete  Properties 

Compressive  strength  f'c 

3,750 

psi 

Modulus  of  elasticity  Ec 

3.7E+6 

psi 

Mass  of  concrete  m0 

4.6584 

lb*sec2/ft 

Unit  weight  yc 

150 

pcf 

Poisson’s  ratio  u 

0.2 

A  typical  design  value  of  compressive  strength  at  28  days  for  a  mass  concrete 
structure  is  2,500  psi.  For  analysis  purpose  it  was  increased  by  a  factor  of  1.5  to 
account  for  aging  effects.  The  modulus  of  elasticity  for  concrete  was 
approximated  by  Equation  3-1  (Ebeling  et  al.  2002): 

Ec=^-54r7a  (3-D 


1  A  table  of  factors  for  converting  non-SI  units  of  measure  to  SI  units  is  found  on  page  xvii. 

2  All  elevations  (el)  cited  in  this  report  are  in  feet  referenced  to  the  National  Geodetic  Vertical 
Datum. 
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where  w  is  the  unit  weight  of  concrete  in  pcf  and  /  'ca  is  the  actual  strength  of 
concrete. 


3.4.  Finite  Element  Models 

Finite  element  models  were  constructed  using  the  preprocessor  incorporated 
within  the  structural  analysis  software  SAP2000  (CS1  1997).  The  generator  bay 
idealized  geometry  was  represented  by  a  hollow  3-D  box  with  five  structural 
elements  labeled  by  numbers  from  1  to  5  as  shown  in  Figure  3-1. 


Upstream 


2 

4 

3 

Base  Slab 

5 

1 

Downstream 


4  ► 


"ground 


a.  Plan  view 


Figure  3-1 .  Idealized  geometry  of  generator  bay  for  composite  finite  element 
model  (Continued) 
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Real  Generator  Bay  Monolith 
Plan  View 


Generator  Bay  Monolith  Idealization 
Plan  View 
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Flow 
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Shear  Wall 
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Equivalent  Inertia  Flexural  Wall 

Flexural  Wall  shear  Wall 


Aground 

Ii,  I2  =  inertia  of  Walls  1  and  2,  respectively 
Ieq  =  equivalent  inertia  of  the  monolith  idealization 
Aeq  =  equivalent  area  of  the  monolith  idealization 


c.  Equivalent  inertia  for  the  upstream  flexural  wall  in  generator  bay  geometry  idealization 


Figure  3-1.  (Concluded) 
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The  idealized  geometry  consists  of  shear  Walls  1  and  2  in  plane  with  the 
upstream/downstream  direction  in  which  the  ground  motion  is  applied.  The 
structural  element  3  is  the  base  slab  of  the  idealized  generator  bay,  and  Walls  4 
and  5  are  the  flexural  walls.  Some  of  the  walls  of  the  water  intake  points  present 
on  the  upstream  side  form  a  channel  through  which  the  intake  water  passes 
(Figure  3-lc).  A  rectangular  flexural  wall  is  an  idealization  of  the  channel  walls 
with  equivalent  inertia  properties  to  those  of  the  real  walls.  Appendix  B  provides 
an  explanation  and  an  example  of  the  calculations  made  for  the  equivalent  area 
and  inertia  for  the  upstream  flexural  wall  for  the  model.  The  dimensions  of  each 
structural  component  of  the  idealized  generator  bay  geometry  are  average  values 
of  the  dimensions  of  the  powerhouses  in  each  group  rounded  to  the  nearest  whole 
number.  The  dimension  for  each  structural  component  of  each  group  is  presented 
in  Table  3-3. 


Table  3-3 

Geometry  of  Composite  Finite  Element  Model 

Group 

Height 

H 

ft 

Upstream/ 

Downstream 

W 

ft 

Out  of  Plane 

B 

ft 

Thickness  f,  ft  | 

Wall  1 

Wall  2 

Base  3 

Wall  4 

Wall  5 

1 

125 

135 

85 

13 

14 

12 

15 

10 

2 

100 

130 

75 

14 

13 

12 

8 

6 

3 

70 

110 

60 

11 

14 

8 

8 

6 

4 

40 

65 

50 

10 

10 

5 

5 

5 

Parametric  studies  were  conducted  to  determine  whether  plane  stress  or  plane 
strain  finite  elements  would  be  used  for  the  analysis.  Plane  stress  finite  elements 
assume  zero  stresses  out  of  plane  of  the  element,  allowing  out-of-plane  deforma¬ 
tion.  The  plane  strain  elements  assume  zero  deformation  out  of  plane  of  the 
element,  permitting  out-of-plane  stress.  Technical  literature  for  concrete  dams 
recommends  plane  stress  elements  because  of  the  construction  joints  between 
monoliths.  This  assumption  permits  deformations  between  monoliths 
(HQUSACE  2003).  Because  powerhouse  (substructure)  monoliths  have  con¬ 
struction  joints  between  them  and  are  constructed  in  a  similar  manner  to  dams, 
plane  stress  elements  were  selected  to  account  for  out-of-plane  deformations  in 
the  structure.  The  result  of  the  parametric  study  between  plane  stress  and  plane 
strain  elements  for  idealized  powerhouse  substructures  is  presented  in 
Appendix  C. 

Nine-node  plane  stress  elements  were  selected  to  represent  the  mass  concrete 
of  the  generator  bay  after  a  parametric  study  was  conducted  of  the  various  types 
of  plane  elements  available  in  SAP2000.  Plane  elements  have  two  degrees  of 
freedom  (vertical  and  horizontal  translation)  at  each  node.  Because  of  the 
formulation  for  the  nine-node  elements,  they  show  better  flexural  behavior  in  the 
model  than  four-node  plane  elements  (Cook  et  al.  2002).  In  the  parametric  study 
between  elements  it  was  found  that  the  fundamental  period  converges  with  fewer 
elements  when  using  nine-node  elements,  reducing  the  number  of  elements 
needed  in  the  model  and  the  time  for  calculations.  The  SAP2000  finite  element 
models  are  shown  in  Figure  3-2  for  the  four  groups.  Appendix  C  contains  an 
explanation  of  the  different  types  of  elements  available  in  SAP2000  and  a 
comparison  of  results  among  them. 
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a.  SAP2000  125-ft-high  composite  model 


b.  SAP2000  100-ft-high  composite  model 


Figure  3-2.  SAP2000  generator  bay  composite  finite  element  models 
(Continued) 
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Figure  3-2.  (Concluded) 
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3.5.  Undamped  Natural  Frequencies  and  Modal 
Mass  Participation  Results 


The  response  to  dynamic  loads  on  the  generator  bay  substructure  will  depend 
on  the  dynamic  properties  of  the  structure  as  well  as  the  frequency  of  the  ground 
motion.  The  natural  periods  and  natural  frequencies  are  two  indices  of  the 
dynamic  structural  properties.  The  undamped  natural  period  is  the  time  required 
for  an  undamped  system  to  complete  one  cycle  of  free  vibration  (Chopra  1995). 
The  natural  frequency  is  inversely  proportional  to  the  natural  period  and 
represents  the  number  of  cycles  in  one  second  (Ebeling  1992). 

A  structure  will  have  a  number  of  vibration  modes,  equal  to  the  number  of 
available  degrees  of  freedom  in  the  finite  element  mesh.  Each  vibration  mode  has 
a  corresponding  period/frequency.  The  total  response  of  the  structure  will  be  a 
combination  of  the  responses  of  each  of  the  modes.  In  a  modal  analysis,  the 
response  of  each  mode  is  calculated  and  combined  to  determine  the  total 
response  of  the  structure  to  dynamic  loads.  The  number  of  modes  needed  for  the 
dynamic  analysis  can  be  determined  from  the  value  of  the  modal  Participating 
Mass  Ratio  (PMR).  The  PMR  of  a  mode  represents  the  contribution  of  the  mode 
to  the  total  response.  For  the  analysis  of  buildings  and  other  structures,  it  is 
recommended  that  the  cumulative  sum  of  the  PMR  be  greater  than  or  equal  to 
0.90  (FEMA  2001).  The  same  criterion  was  used  in  this  study  to  determine  the 
number  of  modes  that  must  be  included  in  the  time-history  modal  analysis. 

A  structural  dynamic  finite  element  analysis  was  performed  on  each  of  the 
four  groups  of  generator  bays  to  calculate  the  first  mode  of  vibration  for  each 
group.  The  analyses  were  done  assuming  a  “dry”  condition  (zero  pool  and 
tailrace  water  elevations)  and  a  “wet”  condition.  The  water  contribution  to  the 
dynamic  analysis  was  represented  by  hydrodynamic  added  mass  on  the  upstream 
and  downstream  sides  of  the  finite  element  models.  A  simplified  procedure, 
Westergaard’s  solution,  was  used  in  the  calculation  of  the  added  water  mass 
(HQUSACE  1995a).  More  information  on  the  calculation  of  hydrodynamic 
added  mass  can  be  found  in  Ebeling  and  Morrison  (1992),  and  HQUSACE 
(1993,  1995b).  Those  calculations  for  the  generator  bays  hydrodynamic  added 
mass  are  presented  in  Appendix  D. 

The  structural  analysis  software  SAP2000  was  used  to  calculate  the  natural 
periods  and  frequencies  of  each  finite  element  model.  The  natural  frequencies 
and  periods  obtained  for  each  group  of  generator  bays  are  presented  in  Table  3-4 
for  the  dry  condition  and  Table  3-5  for  the  wet  condition.  Data  contained  in  these 
tables  show  that  the  models  for  the  wet  condition  exhibited  a  smaller  first  mode 
natural  frequency  than  did  the  models  for  the  dry  condition.  This  observation 
may  be  understood  by  referring  to  the  relationship  for  the  natural  frequency  co 
(rad/sec)  of  a  single  degree  of  freedom  (SDOF)  system  in  Equation  3-2: 


where  k  is  the  stiffness  of  the  structure  and  m  is  the  mass. 
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Table  3-4 

Undamped  Natural  Frequencies  and  Modal  Mass  Participation  Ratio 
for  Generator  Bay  Composite  Models,  Dry  Condition 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

Ti,  sec 

7,  Hz 

Total  No. 

/ 

PMR 

% 

Ti 

sec 

fi 

Hz 

25S9GeD125C5 

125 

0.1008 

9.92 

40 

90.4 

0.007 

143.86 

50S9GeD100C5 

100 

0.0737 

13.58 

34 

90 

0.0069 

143.66 

35S9GeD70C5 

70 

0.048 

20.7 

40 

90 

0.0047 

213.4 

8S9GeD40C5 

40 

0.0286 

35 

26 

91.3 

0.0033 

305.42 

Note:  T,  =  period  for  mode  /;  f  =  frequency  for  mode  1 

;  /  =  mode  number. 

Table  3-5 

Undamped  Natural  Frequencies  and  Modal  Mass  Participation  Ratio 
for  Generator  Bay  Composite  Models,  Wet  Condition 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

Ti,  sec 

7,  Hz 

Total  No. 

/ 

PMR 

% 

T 

sec 

f 

Hz 

25S9GeW125C5 

125 

0.1311 

7.63 

24 

90.5 

0.0127 

78.99 

50S9GeW100C5 

100 

0.095 

10.53 

34 

90.45 

0.0085 

117.32 

35S9GeW70C5 

70 

0.0638 

15.63 

34 

90.7 

0.0069 

145.4 

8S9GeW40C5 

40 

0.0409 

24.46 

16 

90.7 

0.0055 

180.6 

Note:  T,  =  period  for  mode  /;  7  =  frequency  for  mode  1 ;  /  =  mode  number. 

For  a  SDOF  system  with  a  constant  stiffness  k ,  an  increase  in  mass  m  reduces 
co.  Recall  that  the  natural  frequency  /,  in  Flz,  is  given  by  Equation  3-3: 


/ 


CO 

2  n 


(3-3) 


Figure  3-3  presents  the  first  mode  of  vibration  for  the  four  groups  of 
generator  bays  analyzed.  Tables  3-4  and  3-5  summarize  the  first  and  final  modes 
of  vibrations  included  in  each  of  the  analyses  for  the  four  generator  bay  group 
models  as  well  as  the  cumulative  PMR  values  for  each  model.  Figure  3-3  shows 
that  a  shear  type  deformation  behavior,  in  the  first  vibration  mode,  controls  the 
deformation  shape  in  the  models.  This  type  of  behavior  is  found  in  a  structure 
with  small  height  to  width  ratio,  e.g.,  shear  walls. 


3.6.  CEUS  and  PNW  Time-History  Results 

A  dynamic  time-history  analysis  was  performed  on  the  four  groups  of 
generator  bay  substructures  to  determine  the  amplification  of  the  rock  accelera¬ 
tion  due  to  the  ground  motion  at  the  top  of  the  substructure.  A  time-history  of 
relative  acceleration  at  the  top  of  the  substructure  was  calculated  and  then 
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Ext  AM  InM  fn>  rnvlil  nuunnnwlt  _ _ _ 1:1 , jMmgar-^agSSI  BW<  ■«■ cl 

a.  First  mode  of  vibration  for  125-ft-high  composite  model 


b.  First  mode  of  vibration  for  100-ft-high  composite  model 


Figure  3-3.  First  mode  of  vibration  for  generator  bay  composite  models 
(Continued) 
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Figure  3-3.  (Concluded) 


40 
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converted  to  absolute  accelerations  (SA).  To  calculate  the  top-of-substructure 
AFs,  the  SA  at  the  top  of  the  substructure  was  divided  by  the  SA  at  rock  for  each 
frequency.  As  a  result,  a  curve  of  AF  as  a  function  of  frequency  was  created  for 
each  substructure  model  and  for  dry  and  wet  conditions. 

Tables  3-6  and  3-7  list  the  ratio  of  the  total  acceleration  at  a  frequency  of 
50  Flz  (SA50hz)  divided  by  the  PGA  for  dry  and  wet  conditions,  respectively.  The 
PGA  for  the  CEUS  ground  motion  was  0.1 9g  and  0.1 6g  for  the  PNW  ground 
motion.  The  SA  and  AF  results  presented  in  this  chapter  are  for  5  percent  damp¬ 
ing  ratio.  It  is  envisioned  that  these  AF  curves  will  be  used  for  the  decoupled 
analysis  of  powerhouse  superstructures.  The  SA  and  AF  results  for  the  same  top- 
of-substructure  acceleration  time-history  but  processed  for  a  2  percent  damping 
ratio  are  presented  in  Appendix  F.  It  is  envisioned  that  those  AF  curves  will  be 
used  for  the  decoupled  analysis  of  powerhouse  equipment. 


Table  3-6 

Absolute  Acceleration  (SA)  at  Frequency  of  50  Hz  for  Generator 

Bay  Composite  Models,  Dry  Condition 

CEUS 

PNW 

Height 

SA50HZ 

SA/PGA 

SA50HZ 

SA/PGA 

Model  ID 

ft 

g’s 

g’s 

g’s 

g’s 

25S9GeD125C5 

125 

0.61 

3.21 

0.43 

2.69 

50S9GeD100C5 

100 

0.51 

2.68 

0.30 

1.88 

35S9GeD70C5 

70 

0.4 

2.11 

0.23 

1.44 

8S9GeD40C5 

40 

0.46 

2.42 

0.19 

1.19 

Table  3-7 

Absolute  Acceleration  (SA)  at  Frequency  of  50  Hz  for  Generator 

Bay  Composite  Models,  Wet  Condition 

CEUS 

PNW 

Height 

SA50Hz 

SA/PGA 

SA50HZ 

SA/PGA 

Model  ID 

ft 

g’s 

g’s 

g’s 

g’s 

25S9GeW125C5 

125 

0.54 

2.84 

0.41 

2.56 

50S9GeW100C5 

100 

0.56 

2.95 

0.38 

2.38 

35S9GeW70C5 

70 

0.36 

1.89 

0.29 

1.81 

8S9GeW40C5 

40 

0.28 

1.47 

0.20 

1.25 

The  maximum  AFs  at  the  top  of  the  substructure  are  presented  with  their 
corresponding  frequencies  in  Table  3-8  for  the  dry  condition  and  Table  3-9  for 
the  wet  condition.  Note  that  the  maximum  AFs  were  obtained  at  nearly  the  same 
frequency  for  both  ground  motions,  CEUS  and  PNW.  For  the  first  two  groups, 
Group  1  of  125-ft  height  and  Group  2  of  100-ft  height,  the  maximum 
amplification  results  from  the  PNW  ground  motion.  The  substructure  of  the  70-ft 
height  (Group  3)  and  40-ft  height  (Group  4)  has  the  maximum  amplification  in 
response  to  the  CEUS  ground  motion.  The  same  observation  can  be  made  for  wet 
conditions  listed  in  Table  3-9.  For  the  wet  cases,  the  amplifications  were  higher 
than  the  amplification  for  dry  conditions. 
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Table  3-8 

Maximum  AF  and  Its  Corresponding  SA  for  Generator  Bay 
Composite  Models,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

25S9GeD125C5 

125 

9 

2.88 

5.8 

9.9 

2.36 

6.63 

50S9GeD100C5 

100 

13 

2.75 

6.76 

13.6 

1.91 

7.21 

35S9GeD70C5 

70 

20 

1.90 

6.53 

22 

0.79 

3.80 

8S9GeD40C5 

40 

33 

0.75 

3.47 

35 

0.45 

1.18 

Table  3-9 

Maximum  AF  and  Its  Corresponding  SA  for  Generator  Bay 
Composite  Models,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

25S9GeW125C5 

125 

8 

3.27 

7.3 

7.6 

2 

6.04 

50S9GeW100C5 

100 

9 

2.24 

4.84 

10.5 

1.96 

5.85 

35S9GeW70C5 

70 

15 

2.13 

6.34 

15.6 

1.32 

4.79 

8S9GeW40C5 

40 

24.5 

1.16 

4.76 

24.5 

0.48 

2.45 

The  AF  curves  for  each  model,  in  dry  and  wet  conditions,  for  the  two  ground 
motions  are  presented  in  Figures  3-4  to  3-1 1.  Table  3-10  lists  the  figure  number 
for  each  curve.  The  shift  in  maximum  AF  response  to  higher  frequency  (due  to 
the  hydrodynamic  added  mass)  is  demonstrated  in  these  figures.  Note  the  solid 
line  represents  the  dry  condition  SA  and  AF  response  and  the  dashed  line  repre¬ 
sents  the  wet  condition  responses.  The  change  in  response  due  to  the  hydrody¬ 
namic  added  mass  is  dependent  on  the  frequency  change.  Most  of  the  differences 
in  response  between  wet  and  dry  conditions  are  near  the  first  mode  frequency.  As 
the  height  of  structure  decreases,  the  value  for  this  frequency  at  which  maximum 
response  occurs  increases.  An  additional  observation  is  that  the  AF  value  for  wet 
and  dry  conditions  converges  to  nearly  the  same  value  at  higher  frequencies. 


Table  3-10 

Figure  Numbers  of  Generator  Bay  Composite  Model  Results  for 

5  Percent  Damping 

Group 

Height 

ft 

CEUS 

PNW 

1 

125 

Figure  3-4 

Figure  3-8 

2 

100 

Figure  3-5 

Figure  3-9 

3 

70 

Figure  3-6 

Figure  3-10 

4 

40 

Figure  3-7 

Figure  3-11 
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Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


0  5  10  15  20  25  30  35  40  45  50 


Frequency  (Hz) 

b.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


Frequency  (Hz) 

c.  CEUS  ground  motion  response  spectrum 

Figure  3-4.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
125-ft  height  (Dry  Model  ID:  25S9GeD125C5,  Wet  Model  ID: 
25S9GeW125C5),  for  CEUS  ground  motion  with  5  percent  damping 
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Figure  3-5.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
100-ft  height  (Dry  Model  ID:  50S9GeD100C5,  Wet  Model  ID: 
50S9GeW100C5),  for  CEUS  ground  motion  with  5  percent  damping 
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Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


b.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


Figure  3-6.  Response  spectra  and  AF  curves,  generator  bay,  composite  model 
70-ft  height  (Dry  Model  ID:  35S9GeD70C5,  Wet  Model  ID: 
35S9GeW70C5),  for  CEUS  ground  motion  with  5  percent  damping 


Chapter  3  Time-History  Finite  Element  Results  for  the  Generator  Bay  Monoliths 


SA  (g)  q-  SA  (g)  (Amplification  Factor) 


Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 
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Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 
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Frequency  (Hz) 

c.  CEUS  ground  motion  response  spectrum 

Figure  3-7.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
40-ft  height  (Dry  Model  ID:  8S9GeD40C5,  Wet  Model  ID: 
8S9GeW40C5),  for  CEUS  ground  motion  with  5  percent  damping 
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Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


Frequency  (Hz) 

b.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


Frequency  (Hz) 

c.  PNW  ground  motion  response  spectrum 

Figure  3-8.  Response  spectra  and  AF  curves,  generator  bay,  composite  model 
125-ft  height  (Dry  Model  ID:  25S9GeD125C5,  Wet  Model  ID: 
25S9GeW125C5),  for  PNW  ground  motion  with  5  percent  damping 
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a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 
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.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 
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c.  PNW  ground  motion  response  spectrum 

Figure  3-9.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
100-ft  height  (Dry  Model  ID:  50S9GeD100C5,  Wet  Model  ID: 
50S9GeW100C5),  for  PNW  ground  motion  with  5  percent  damping 
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Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 
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b.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


Frequency (Hz) 

c.  PNW  ground  motion  response  spectrum 

Figure  3-10.  Response  spectra  and  AF  curves,  generator  bay,  composite  model 
70-ft  height  (Dry  Model  ID:  35S9GeD70C5,  Wet  Model  ID: 
35S9GeW70C5),  for  PNW  ground  motion  with  5  percent  damping 
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Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


Frequency  (Hz) 

.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


Frequency (Hz) 

c.  PNW  ground  motion  response  spectrum 

Figure  3-11.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
40-ft  height  (Dry  Model  ID:  8S9GeD40C5,  Wet  Model  ID: 
8S9GeW40C5),  for  PNW  ground  motion  with  5  percent  damping 
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3.7.  Response  Amplification  Sensitivity  to 
Location  along  the  Top  of  the  Substructure 
for  the  Generator  Bay  Composite  Model 


The  results  presented  in  Section  3.5  were  obtained  by  calculating  the  spectral 
acceleration  at  the  top  midnode  of  the  composite  finite  element  model.  A  para¬ 
metric  study  was  performed  to  investigate  the  variation  in  amplification  across 
the  top  of  the  substructure.  The  spectral  accelerations  (5  percent  damping)  and 
AFs  versus  frequency  curves  were  calculated  at  the  three  nodes  identified  in 
Figure  3-12.  The  finite  element  model  used  for  this  analysis  was  the  125-ft- 
height  composite  model,  wet  condition. 


Figure  3-12.  Location  of  the  three  nodes  in  the  125-ft-height  generator  bay  composite  model  used  to 
calculate  AFs 

Table  3-11  lists  values  of  SA  at  50  FIz  and  its  ratio  to  the  PGA  for  the  two 
ground  motions  of  CEUS  and  PNW.  Table  3-12  lists  the  maximum  AF  and  the 
corresponding  frequency  at  the  three  different  locations  along  the  top  of  the 
substructure.  Figure  3-13  shows  the  variation  in  spectral  acceleration  with  fre¬ 
quency  for  the  three  locations  at  top  of  the  substructure,  and  Figure  3-14  shows 
the  corresponding  AF  curves.  Because  of  the  influence  of  the  added  water  mass, 
as  explained  in  Appendix  D,  the  responses  near  the  upstream  side  are  higher  than 
those  for  the  midnode  and  the  downstream  node.  At  frequencies  near  the  first 
mode  of  vibration  of  the  finite  element  model,  the  AF  at  the  middle  of  the 
substructure  is  equal  to  the  AF  at  the  downstream  side. 
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Table  3-11 

SA  at  50  Hz  at  Three  Different  Locations  along  the  Top  of  the 
125-ft-Height  Composite  Substructure  Model  for  the  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g's 

SA/ PGA 
g's 

SA50HZ 

g's 

SA/ PGA 
g's 

25S9GeW125C5- 

Upstream 

125 

0.73 

3.84 

0.5 

3.13 

25S9GeW125C5- 

Midnode 

125 

0.54 

2.84 

0.41 

2.56 

25S9GeW125C5- 

Downstream 

125 

0.56 

2.95 

0.41 

2.56 

Table  3-12 

Maximum  AF  and  Its  Corresponding  Frequency  at  Three  Different 
Locations  along  the  Top  of  the  125-ft-Height  Composite 
Substructure  Model  for  the  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g's 

(AF)sa 
at  f 

f 

Hz 

SA 

g's 

(AF)sA 
at  f 

25S9GeW125C5- 

Upstream 

125 

7.62 

3.55 

8.45 

7.62 

2.32 

7.09 

25S9GeW125C5- 

Midnode 

125 

8 

3.27 

7.3 

7.6 

2 

6.04 

25S9GeW125C5- 

Downstream 

125 

7.62 

3.08 

7.33 

7.62 

2.01 

6.15 

3.8.  Response  Amplification  Sensitivity  to  Shear 
Wall  Thickness  for  the  Generator  Bay  Composite 
Finite  Element  Models 

The  influence  of  shear  wall  thickness  of  Wall  1  plus  Wall  2  on  the  computed 
results  for  the  generator  bay  composite  model  was  investigated.  The  125-ft- 
height  model  in  the  wet  condition  was  used  for  the  additional  pair  of  analyses. 

The  total  thickness  of  Walls  1  plus  2  in  the  original  finite  element  model  is  27  ft. 
For  this  parametric  study,  the  shear  wall  thickness  of  the  standard  finite  element 
model  was  decreased  and  increased  by  50  percent,  resulting  in  a  total  thickness  of 
13.5  ft  and  40.5  ft,  respectively. 

Table  3-13  lists  the  first  mode  period  and  frequency,  the  number  of  modes 
used  in  the  analyses,  and  the  PMR.  Table  3-14  lists  the  SA  at  50  Hz  and  its  ratio 
to  the  PGA  for  the  two  ground  motions.  Table  3-15  lists  the  maximum  AFs  for 
the  three  shear  wall  thickness  considered. 

The  results  show  a  change  in  the  model  stiffness  due  to  change  in  shear  wall 
thickness.  The  fundamental  mode  frequency  increases  by  1.11  Hz  when  the  total 
shear  wall  thickness  is  increased  from  13.5  ft  to  27  ft.  Conversely,  the  fundamen¬ 
tal  mode  frequency  decreases  by  0.52  Hz  when  the  total  shear  wall  thickness  is 
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Frequency  (Hz) 


- Upstream 

.  Midnode 

— * —  Downstream 


a.  CEUS  Ground  Motion 


Frequency  (Hz) 


- Upstream 

.  Midnode 

— * —  Downstream 


b.  PNW  Ground  Motion 


Figure  3-13.  Response  spectrum  curves  for  nodes  at  three  different  locations 
along  the  top-of-substructure  of  the  generator  bay  125-ft-height, 
composite  model,  in  the  wet  condition,  using  5  percent  damping 
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(Amplification  Factor)SA  (Amplification  Factor)SA 


Figure  3-14.  AF  curves  for  nodes  at  three  different  locations  along  the  top-of- 

substructure  of  the  generator  bay,  125-ft-height,  composite  model  in 
the  wet  condition,  using  5  percent  damping 
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Table  3-13 

First  Mode  of  Vibration  and  PMR  for  Three  125-ft-Height  Generator 
Bay  Composite  Models 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

sec 

u 

Hz 

Total  No. 

/ 

PMR 

% 

Ti 

sec 

f, 

Hz 

25S9GeW125C5-13.5  ft 

125 

0.1535 

6.52 

24 

90 

0.0135 

73.89 

25S9GeW125C5-27  ft 

125 

0.1311 

7.63 

24 

90.5 

0.0127 

78.99 

25S9GeW125C5-40.5  ft 

125 

0.1227 

8.15 

24 

91.4 

0.0119 

84.16 

Table  3-14 

SA  at  50  Hz  for  the  Three  Finite  Element  Models  Using  the  CEUS 
and  PNW  Ground  Motions 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50hz 

g's 

SA / PGA 
g's 

SA50hz 

g's 

SA / PGA 
g's 

25S9GeW125C5-13.5  ft 

125 

0.49 

2.58 

0.4 

2.50 

25S9GeW125C5-27  ft 

125 

0.54 

2.84 

0.41 

2.56 

25S9GeW125C5-40.5  ft 

125 

0.6 

3.16 

0.40 

2.50 

Table  3-15 

Maximum  AF  and  Its  Corresponding  Frequency  for  the  Three  125-ft- 
Height  Generator  Bay  Models 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g's 

(AF)sa 
at  f 

f 

Hz 

SA 

g's 

(AF)sa 
at  f 

25S9GeW125C5-13.5  ft 

125 

6.52 

2.84 

7.14 

6 

2.03 

5.78 

25S9GeW125C5-27  ft 

125 

7.5 

3.01 

7.1 

7.62 

2 

6.12 

25S9GeW125C5-40.5  ft 

125 

8.15 

3.80 

8.26 

8.15 

1.63 

5.14 

decreased  from  40.5  ft  to  27  ft.  Figure  3-15  shows  the  absolute  acceleration 
curves  for  the  three  models,  and  Figure  3-16  shows  the  AF  curves  for  the  two 
ground  motions  considered.  The  computed  results  show  that  the  shear  wall 
thickness  can  have  a  significant  influence  on  the  computed  results,  depending 
upon  the  frequency  content  of  the  ground  motion. 
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IB)  VS  (6)  VS 


Figure  3-15.  Response  spectrum  curves  at  top-of-substructure  for  three  different 
shear  wall  thicknesses  (sum  of  Wall  1  and  2  thicknesses)  of  the 
generator  bay,  125-ft-height,  composite  model,  in  the  wet  condition, 
using  5  percent  damping 
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(Amplification  Factor)SA  (Amplification  Factor)SA 


Frequency  (Hz) 


a.  CEUS  Ground  Motion 


Frequency  (Hz) 
b.  PNW  Ground  Motion 

Figure  3-16.  AF  curves  at  top-of-substructure  for  three  different  shear  wall 

thicknesses  (sum  of  Wall  1  and  2  thicknesses)  of  the  generator  bay, 
125-ft-height,  composite  model,  in  the  wet  condition,  using  5  percent 
damping 
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4  Time-History  Finite 
Element  Results  for  the 
Erection/Service  Bay 
Monolith 


4.1.  Introduction 

This  chapter  discusses  the  results  from  the  analysis  of  a  series  of  time- 
history/modal  response  of  two  models  of  erection/service  bay  substructure  mono¬ 
liths  to  earthquake  loadings.  Two  acceleration  time-histories  were  used  in  this 
study,  one  representing  a  CEUS  earthquake  event  and  a  second  representing  a 
PNW  event.  A  block  model  and  a  block-frame-shear  wall  finite  element  model  of 
each  idealized  powerhouse  substructure  were  created  using  the  structural  analysis 
software  SAP2000  (CSI  1997).  The  acceleration  time-history  was  specified  at  the 
base  of  the  models,  and  time-histories  at  the  top  of  the  substructure  were  calcu¬ 
lated  and  then  converted  into  response  spectra.  The  response  spectra  at  the  top  of 
the  substructure  were  converted  to  AFs  versus  frequency  curves  for  use  in  a 
decoupled  response  spectra  modal  analysis.  The  AF  curves  presented  in  this 
chapter  are  plotted  up  to  50  Hz.  However,  the  ground  motions  used  in  this  study 
have  frequency  content  up  to  33  Hz.  The  use  of  AF  curves  above  33  Hz  is  not 
recommended. 


4.2.  Geometry  of  Four  Height  Monoliths 

Fifty-four  cross-sections  of  the  erection/service  bay  substructures  of  Corps 
powerhouses  were  reviewed.  The  erection/service  bay  substructures  were 
assigned  to  one  of  four  groups  according  to  their  heights.  Average  heights, 
upstream/downstream  dimensions,  plan  widths,  pool  water  elevation,  and  tailrace 
elevation  were  determined  for  each  substructure  group.  Appendix  A  summarizes 
the  Corps  erection/service  bay  monoliths  geometry  reviewed  in  this  study. 

Two  types  of  finite  element  models  were  used  to  idealize  the  geometry  of 
the  erection/service  bay  substructure.  One  is  referred  to  as  a  block  model  (Fig¬ 
ure  4-1),  and  the  second  as  a  block-frame-shear  wall  model  (Figure  4-2).  The 
block  idealization  represents  a  stiff  model  of  the  substructure,  and  the  block- 
frame-shear  wall  model  represents  a  more  flexible  substructure. 
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Figure  4-1 .  Geometry  of  erection/service  bay  block  finite  element  model 


Figure  4-2.  Geometry  of  an  erection/service  bay  block-frame-shear  wall  finite  element  model 
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The  dimensions  for  each  model  are  presented  in  Table  4-1.  The  shear  wall  thick¬ 
ness  for  the  block- frame-shear  wall  model  was  selected  by  inspection  of  the 
construction  drawings  for  different  Corps  powerhouses.  An  average  8-ft  wall 
thickness  was  judged  to  represent  the  thickness  of  the  shear  wall(s). 


Table  4-1 

Geometry  of  Block  and  Block-Frame-Shear  Wall  Finite  Element 
Models 

Group 

Height 

H 

ft 

Upstream/ 

downstream 

W 

ft 

Out  of  Plane 
Width 

B 

ft 

Pool  Water 

El 

ft 

Tailrace 

El 

ft 

Range  in 
Height  for 
Group 
ft 

1 

110 

125 

75 

130 

75.5 

90-135 

2 

75 

120 

60 

95 

49 

60-90 

3 

45 

110 

65 

65 

35 

30-55 

4 

20 

60 

40 

40 

16 

0-30 

4.3.  Material  Properties 

The  material  used  in  the  construction  of  powerhouse  substructures  is  con¬ 
crete.  For  mass  concrete  2,500  psi  of  compressive  strength  is  typically  selected  as 
the  design  value  and  3,000  psi  for  the  reinforced  concrete  frames.  To  account  for 
aging  effects,  the  compressive  strength  for  both  entities  was  increased  by  a  factor 
of  1.5.  The  general  properties  for  mass  concrete  are  listed  in  Table  4-2  and  the 
concrete  properties  for  the  frame  elements  given  in  Table  4-3.  The  value  assigned 
to  Young’s  modulus  of  the  mass  and  structural  concrete  was  computed  using 
Equation  3-1  for  the  aged  compressive  strength  f'  values. 


Table  4-2 

Material  Properties  Used  in  the  Powerhouse  Plane-Stress 

Continua  Finite  Elements 

Mass  Concrete  Properties 

Compressive  strength  f  'c 

3,750 

psi 

Modulus  of  elasticity  Ec 

3.7E+6 

psi 

Mass  of  concrete  mc 

4.6584 

lb*sec2/ft 

Unit  weight  yc 

150 

pcf 

Poisson’s  ratio  o 

0.2 

Table  4-3 

Material  Properties  Used  in  the  Frame  Elements 

Concrete  Properties 

Compressive  strength  f  'c 

4,500  psi 

Modulus  of  elasticity  Ec 

4.1  E+6  psi 

Mass  of  concrete  m0 

4.6584  lb*sec2/ft 

Unit  weight  yc 

150  pcf 

Poisson’s  ratio  u 

0.2 
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4.4.  Finite  Elements  Models 


Nine-node  plane  stress  elements  were  selected  to  represent  the  mass  of 
concrete  in  the  substructure  for  the  block  and  block-frame-shear  wall  models. 

The  plane  stress  elements  assume  zero  stresses  out  of  plane  to  the  elements  but 
allow  for  out-of-plane  deformation.  This  assumption  is  consistent  with  the 
recommendation  in  the  technical  literature  (HQUSACE  2003)  because  of  the 
presence  of  construction  joints  between  monoliths. 

In  the  block- frame-shear  wall  model,  flexural  deformation  is  expected 
because  of  the  frame’s  rotational  degrees  of  freedom  inside  the  model.  The  use  of 
nine-node  elements  reduces  the  shear  locking  behavior  of  lower  order  elements 
(Cook  et  al.  2002).  Shear  locking  is  a  behavior  that  allows  an  element  to  over¬ 
estimate  its  shear  deformations.  When  this  occurs,  more  forces  are  needed  to 
deform  the  structure  in  flexure  (i.e.,  an  analytically  stiff  element). 

To  represent  the  internal  frame  in  the  substructure,  frame  elements  with  three 
degrees  of  freedom  (two  translations  and  one  rotation)  at  each  node  were 
selected.  The  frame  elements  were  extended  into  the  plane  finite  elements  of  the 
mass  concrete  to  avoid  the  development  of  analytical  hinges  at  their  connection. 

SAP2000  structural  analysis  software  was  used  to  create  the  models  for  the 
erection/service  bay  substructure.  Figure  4-3  shows  the  four  groups  of  block 
models.  The  minimum  number  of  elements  was  determined  by  increasing  the 
refinement  of  the  mesh  until  a  constant  first  mode  period  value  was  computed. 
Figure  4-4  shows  the  SAP2000  models  for  the  block- frame-shear  wall.  To 
account  for  cracking  effects  in  the  frame  elements,  the  gross  inertia  on  these 
elements  was  reduced  by  20  percent  (Ebeling  et  al.  2002). 


4.5.  Undamped  Natural  Frequencies  and  Modal 
Mass  Participation  Results 

Undamped  natural  frequencies  and  periods  for  the  block  models  are  reported 
in  Table  4-4  for  the  dry  condition  and  Table  4-5  for  the  wet  condition.  A  5  per¬ 
cent  damping  was  assumed  for  the  substructure  in  this  chapter.  The  AF  curves 
presented  in  this  chapter  were  obtained  using  response  spectra  curves  for  5  per¬ 
cent  damping.  The  SA  and  AF  results  for  2  percent  damping  are  presented  in 
Appendix  F  for  use  with  equipment.  The  number  of  modes  considered  for  each 
model  was  increased  until  a  cumulative  PMR  over  0.90  was  attained.  Tables  4-4 
and  4-5  also  list  the  period  and  frequency  for  the  final  mode  considered  for  each 
substructure  model.  The  modal  analysis  results  showed  that  the  hydrodynamic 
added  mass  reduced  the  natural  frequency  of  the  model,  making  the  model 
behave  as  if  it  is  more  “flexible.”  This  effect  is  more  evident  in  the  substructure 
groups  of  1 1 0-ft-height  and  75-ft-height  substructure  where  the  difference  in  SA 
and  AF  peak  responses  between  dry  and  wet  models  is  greater  than  for  sub¬ 
structure  groups  3  and  4. 
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a.  110-ft  height 


b.  75-ft  height 


Figure  4-3.  Erection/service  bay  block  finite  element  models  (Continued) 
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Figure  4-3.  (Concluded) 
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a.  110-ft  height 


b.  75-ft  height 


Figure  4-4.  Erection/service  bay  block-frame-shear  wall  finite  element  models  (Continued) 
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Figure  4-4.  (Concluded) 
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Table  4-4 

Undamped  Natural  Frequencies  and  Modal  Mass  Participation 

Ratio  for  Erection/Service  Bay  Block  Models,  Dry  Condition 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

71,  sec 

fu  Hz 

Total  No. 

/ 

PMR 

% 

T, 

sec 

f, 

Hz 

16S9ErD110B5 

110 

0.0915 

10.93 

8 

92.3 

0.0155 

64.45 

16S9ErD75B5 

75 

0.0554 

18.04 

8 

90.2 

0.0125 

79.69 

16S9ErD45B5 

45 

0.0304 

32.91 

12 

90.6 

0.0069 

144.95 

8S9ErD20B5 

20 

0.0131 

76.13 

12 

90.4 

0.0036 

274.19 

Note:  T,  =  period  for  mode  /;  f,  =  frequency  for  mode  1 ;  /  =  mode  number. 

Table  4-5 

Undamped  Natural  Frequencies  and  Modal  Mass  Participation 

Ratio  for  Erection/Service  Bay  Block  Models,  Wet  Condition 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode  j 

71,  sec 

fu  Hz 

Total  No. 

/ 

PMR 

% 

Ti 

sec 

f, 

Hz 

16S9ErW110B5 

110 

0.1038 

9.63 

9 

92.5 

0.0189 

52.75 

16S9ErW75B5 

75 

0.0627 

15.96 

10 

90.7 

0.0132 

75.78 

16S9ErW45B5 

45 

0.0357 

27.99 

14 

90.8 

0.0077 

129.9 

8S9ErW20B5 

20 

0.0213 

46.8 

14 

90.1 

0.004 

250.77 

Note:  71  =  period  for  mode  i;  7  =  frequency  for  mode  1 ;  /  =  mode  number. 

The  first  mode  of  vibration  of  the  four  groups  of  block  models  is  shown  in 
Figure  4-5.  A  shear  deformation  behavior  is  dominant  for  these  models  as 
expected  for  a  structure  with  a  low  height  to  width  ratio  subjected  to  lateral  loads. 

Tables  4-6  and  4-7  list  the  periods  and  frequencies  for  the  block- frame-shear 
wall  models  for  dry  and  wet  conditions,  respectively,  and  their  corresponding 
accumulated  participating  mass  ratio.  The  same  hydrodynamic  added  mass  effect 
on  the  natural  periods  is  observed  in  the  erection/service  bay  substructure  models 
as  for  the  generator  bay,  explained  in  Section  3.4  of  this  report. 

The  block- frame-shear  wall  models  showed  more  flexibility  than  the  block 
models  because  of  the  internal  frame  in  the  substructure.  Figure  4-6  shows  the 
first  mode  of  vibration  of  the  four  groups  of  the  block- frame-shear  wall  models. 


4.6.  CEUS  and  PNW  Time-History  Results  for 
Block  Finite  Element  Models 

4.6.1  Analysis  results 

A  time-history  analysis  was  performed  on  the  four  erection/service  bay  sub¬ 
structure  block  models.  The  amplified  acceleration  time-histories  were  recorded 
at  the  midtop  node  of  the  finite  element  mesh.  These  time-histories  were  then 
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b.  75-ft  height 


Figure  4-5.  First  mode  of  vibration  for  erection/service  bay  block  models  (Continued) 
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c.  45-ft  height 


d.  20-ft  height 


Figure  4-5.  (Concluded) 
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Table  4-6 

Undamped  Natural  Frequencies  and  Modal  Mass  Participation 

Ratio  for  Erection/Service  Bay  Block-Frame-Shear  Wall  Models, 

Dry  Condition 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

7"i,  sec 

f,,  Hz 

Total  No. 

/ 

PMR 

% 

Ti 

sec 

f, 

Hz 

14S9ErD110FW5 

110 

0.0964 

10.37 

14 

90.39 

0.0136 

73.72 

20S9ErD75FW5 

75 

0.0624 

16.02 

15 

90.06 

0.0098 

101.8 

12S9ErD45FW5 

45 

0.0362 

27.04 

15 

90.35 

0.0075 

132.94 

8S9ErD20FW5 

20 

0.0135 

74.05 

20 

91.33 

0.0025 

405.21 

Table  4-7 

Undamped  Natural  Frequencies  and  Modal  Mass  Participation 

Ratio  for  Erection/Service  Bay  Block-Frame-Shear  Wall  Models, 

Wet  Condition 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

Ti,  sec 

fr  Hz 

Total  No. 

/ 

PMR 

% 

Ti 

sec 

f, 

Hz 

14S9ErW110FW5 

110 

0.1193 

8.38 

18 

93.3 

0.0141 

71.03 

14S9ErW75FW5 

75 

0.0773 

12.94 

15 

90.54 

0.0112 

89.58 

12S9ErW45FW5 

45 

0.0479 

20.89 

20 

92.76 

0.0071 

141.35 

8S9ErW20FW5 

20 

0.0281 

35.52 

30 

91.88 

0.0027 

372.33 

converted  to  SA  at  top-of-substructure  for  5  percent  damping.  To  develop  the  AF 
curves,  the  SA  at  top-of-substructure  was  divided  by  the  SA  for  5  percent  damp¬ 
ing  at  rock  for  the  same  frequency. 

Tables  4-8  and  4-9  list  the  SA  at  top-of-substructures  for  dry  and  wet  condi¬ 
tions,  respectively,  at  a  frequency  of  50  Flz  and  the  ratio  with  the  PGA  for  each 
of  the  two  ground  motions.  The  in-rock  PGA  for  the  CEUS  ground  motion  was 
taken  as  0. 19g  and  for  the  PNW  ground  motion,  0. 16g.  The  data  given  in  these 
two  tables  show  that  the  AF  values  decrease  with  the  height  of  substructure. 

The  maximum  AFs  and  their  corresponding  frequency  are  listed  in 
Table  4-10  for  the  dry  condition  and  Table  4-11  for  the  wet  condition  for  the  two 
ground  motions.  The  results  show  that  the  20-ft-height  block  models  can  be 
considered  as  rigid  structures  because  the  AF  is  near  one  for  dry  and  wet  condi¬ 
tions  for  both  earthquakes.  This  was  expected  because  its  natural  frequency  is 
well  beyond  the  maximum  ground  motion  frequency  contents  of  33  Flz  (as 
reflected  in  the  zero  Fourier  amplitudes  in  Figures  2-9  and  2-10). 

Figures  4-7  to  4-14  show  the  AF  curves,  the  response  spectrum  at  top-of- 
substructure,  and  the  response  spectrum  at  rock  for  the  block  models  with 
5  percent  damping.  The  results  are  presented  first  for  the  CEUS  ground  motion 
and  then  for  the  PNW  ground  motion.  Table  4-12  lists  the  figure  numbers  of  SA 
and  AF  results  on  this  section. 
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E3  Mode  1  Period  0.0964  seconds 


a.  110-ft  height 


Click  on  Joint  for  modal  components 

Start  Animation  |  <J=>  <=t>  |lb-ft 

▼ 

b.  75-ft  height 

Figure  4-6.  First  mode  of  vibration  for  erection/service  bay  block-frame-shear  wall  models  (Continued) 
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Figure  4-6.  (Concluded) 
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Table  4-8 

SA  at  Frequency  of  50  Hz  for  Erection/Service  Bay  Block  Models, 

Dry  Condition 

CEUS 

PNW 

Height 

SA50HZ 

SA/PGA 

SA50HZ 

SA/PGA 

Model  ID 

ft 

g’s 

g’s 

g’s 

g’s 

16S9ErD110B5 

110 

0.58 

3.05 

0.41 

2.56 

16S9ErD75B5 

75 

0.4 

2.11 

0.22 

1.38 

16S9ErD45B5 

45 

0.27 

1.42 

0.17 

1.06 

8S9ErD20B5 

20 

0.2 

1.05 

0.16 

1 

Table  4-9 

SA  at  Frequency  of  50  Hz  for  Erection/Service  Bay  Block  Models, 

Wet  Condition 

CEUS 

PNW 

Height 

SA50HZ 

SA/PGA 

SA50HZ 

SA/PGA 

Model  ID 

ft 

g’s 

g’s 

g’s 

g’s 

16S9ErW110B5 

110 

0.68 

3.58 

0.45 

2.81 

16S9ErW75B5 

75 

0.39 

2.05 

0.31 

1.94 

16S9ErW45B5 

45 

0.27 

1.42 

0.18 

1.13 

8S9ErW20B5 

20 

0.19 

1.00 

0.16 

1 

Table  4-10 

Maximum  Amplification  Factor  (AF)Sa  and  Its  Corresponding  SA  for 
Erection/Service  Bay  Block  Models,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g’s 

(AF)sA 

at  f 

f 

Hz 

SA 

g’s 

(AF)sA 
at  f 

16S9ErD110B5 

110 

11 

2.43 

5.76 

11 

1.98 

5.97 

16S9ErD75B5 

75 

18.04 

2.08 

6.16 

16.5 

0.99 

4.38 

16S9ErD45B5 

45 

32.91 

1.07 

4.95 

4 

0.45 

1.11 

8S9ErD20B5 

20 

- 

- 

- 

- 

- 

- 

Table  4-11 

Maximum  Amplification  Factor  (AF)Sa  and  Its  Corresponding  SA  for 
Erection/Service  Bay  Block  Models,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

16S9ErW110B5 

110 

9 

3.35 

6.76 

9.63 

2.45 

6.95 

16S9ErW75B5 

75 

15 

2.11 

6.26 

15.96 

1.44 

5.11 

16S9ErW45B5 

45 

28 

0.88 

3.55 

4 

0.45 

1.47 

8S9ErW20B5 

20 

- 

- 

- 

- 

- 

- 
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Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 
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c.  CEUS  ground  motion  response  spectrum 

Figure  4-7.  Response  spectra  and  AF  curves,  erection/service  bay,  block  model 
1 1 0-ft  height  (Dry  Model  ID:  1 6S9ErD1 1 0B5,  Wet  Model  ID: 

1 6S9ErW1 1 0B5),  for  CEUS  ground  motion  with  5  percent  damping 
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c.  CEUS  ground  motion  response  spectrum 

Figure  4-8.  Response  spectra  and  AF  curves,  erection/service  bay,  block  model, 
75-ft  height  (Dry  Model  ID:  16S9ErD75B5,  Wet  Model  ID: 
16S9ErW75B5),  for  CEUS  ground  motion  with  5  percent  damping 
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c.  CEUS  ground  motion  response  spectrum 

Figure  4-9.  Response  spectra  and  AF  curves,  erection/service  bay,  block  model 
45-ft  height  (Dry  Model  ID:  16S9ErD45B5,  Wet  Model  ID: 
16S9ErW45B5),  for  CEUS  ground  motion  with  5  percent  damping 
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Figure  4-10.  Response  spectra  and  AF  curves,  erection/service  bay,  block 
model,  20-ft  height  (Dry  Model  ID:  8S9ErD20B5,  Wet  Model  ID: 
8S9ErW20B5),  for  CEUS  ground  motion  with  5  percent  damping 
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Figure  4-11.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  110-ft  height  (Dry  Model  ID:  16S9ErD110B5,  Wet  Model  ID: 
16S9ErW1 10B5),  for  PNW  ground  motion  with  5  percent  damping 
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Figure  4-12.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  75-ft  height  (Dry  Model  ID:  16S9ErD75B5,  Wet  Model  ID: 
16S9ErW75B5),  for  PNW  ground  motion  with  5  percent  damping 
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Figure  4-13.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  45-ft  height  (Dry  Model  ID:  16S9ErD45B5,  Wet  Model  ID: 
16S9ErW45B5),  for  PNW  ground  motion  with  5  percent  damping 
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Figure  4-14.  Response  spectra  and  AF  curves,  erection/service  bay,  block 
model,  20-ft  height  (Dry  Model  ID:  8S9ErD20B5,  Wet  Model  ID: 
8S9ErW20B5),  for  PNW  ground  motion  with  5  percent  damping 
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Table  4-12 

Figure  Numbers  of  Erection/Service  Bay  Block  Model  Results  for 

5  percent  damping 

Group 

Height 

CEUS 

PNW 

1 

110 

Figure  4-7 

Figure  4-11 

2 

75 

Figure  4-8 

Figure  4-12 

3 

45 

Figure  4-9 

Figure  4-13 

4 

20 

Figure  4-10 

Figure  4-14 

4.6.2  Response  amplification  sensitivity  to  location  at  top  of  the 
substructure  for  the  erection/service  bay  block  model 

To  investigate  the  variation  of  the  results  due  to  the  location  at  top  of  the 
substructure  of  the  erection  bay  block  model,  the  absolute  accelerations 
(5  percent  damping)  and  AFs  versus  frequency  curves  were  calculated  at  the 
three  nodes  identified  in  Figure  4-15.  The  finite  element  model  used  for  this 
analysis  was  the  1 1 0-ft-height  erection/service  bay  block  model,  in  wet 
condition. 


Figure  4-15.  Location  of  the  three  nodes  in  the  1 1 0-ft-height  erection/service  bay  block  model  used  to 
calculate  the  AFs 
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Table  4-13  lists  the  SA  at  50  Hz  and  its  ratio  with  the  PGA  at  three  locations 
in  the  substructure.  Table  4-14  lists  the  maximum  AF  and  its  frequency  at  three 
locations  in  the  substructure.  Figure  4-16  shows  the  SA  for  the  three  locations  at 
top  of  the  substructure  for  both  the  CEUS  and  the  PNW  ground  motions. 

Figure  4-17  presents  the  AF  curves  at  the  three  different  locations  along  the  top 
of  the  substructure. 


Table  4-13 

SA  at  50  Hz  at  Three  Different  Locations  along  the  Top  of  the 
110-ft-Height  Block  Substructure  Model  for  the  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g's 

SA5ohz/PGA 

g's 

SA50HZ 

g's 

SA5ohz/PGA 

g's 

1 6S9ErW1 1 0B5-Upstream 

110 

0.79 

4.16 

0.48 

3.00 

16S9ErW1 1 0B5-Midnode 

110 

0.68 

3.58 

0.45 

2.81 

1 6S9ErW1 1 0B5-Downstream 

110 

0.68 

3.58 

0.45 

2.81 

Table  4-14 

Maximum  AF  and  Its  Corresponding  Frequency  at  Three  Different 
Locations  along  the  Top  of  the  110-ft-Height  Block  Substructure 

Model  for  the  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g's 

(AF)sA 
at  f 

f 

Hz 

SA 

g's 

(AF)sA 
at  f 

1 6S9ErW1 1 0B5-Upstream 

110 

9 

3.64 

7.37 

9.62 

2.72 

7.8 

16S9ErW1 1 0B5-Midnode 

110 

9 

3.35 

6.76 

9.63 

2.45 

6.95 

1 6S9ErW  1 1 0B5-Downstream 

110 

9 

3.43 

6.94 

9.62 

2.55 

7.32 

As  shown  in  Figures  4-16  and  4-17  for  frequencies  near  the  first  mode  of 
vibration  for  the  block  model,  the  SA  and  the  AF  are  similar  at  the  midnode  and 
the  downstream  node.  The  values  at  the  upstream  side  are  slightly  higher  because 
of  the  influence  of  the  additional  water  mass  included  in  the  upstream  top  node  as 
explained  in  Appendix  D. 


4.7.  CEUS  and  PNW  Time-History  Results  for 
Block-Frame-Shear  Wall  Finite  Element  Models 

4.7.1  Analysis  Results 

To  account  for  the  internal  frames  in  the  erection/service  bay  substructure,  a 
block- frame-shear  wall  finite  element  model  was  created  and  analyzed  using  the 
two  ground  motions  of  this  study.  It  is  observed  in  computed  results  that  the 
internal  frame  makes  the  model  more  flexible  than  the  block  models,  increasing 
the  amplification  at  top  of  the  substructure. 
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Figure  4-16.  Response  spectrum  curves  for  nodes  at  three  different  locations 

along  the  top-of-substructure  of  the  erection  bay  110-ft-height  block 
model,  in  wet  condition,  using  5  percent  damping 
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Figure  4-17.  AF  curves  for  nodes  at  three  different  locations  along  the  top-of- 
substructure  of  the  erection  bay,  110-ft-height  block  model,  in  wet 
condition,  using  5  percent  damping 
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Tables  4-15  and  4-16  list  the  SA  at  50  Hz  for  the  two  ground  motions  and  its 
ratio  to  the  in-rock  PGA  for  the  wet  and  dry  conditions,  respectively.  The  results 
in  this  section  were  obtained  assuming  5  percent  damping  for  the  substructure 
model.  The  AF  curves  were  obtained  using  response  spectrum  curves  for 
5  percent  damping  and  used  in  the  decoupled  analysis  of  the  superstructure.  The 
results  for  2  percent  damping  are  presented  in  Appendix  F  for  use  with  the 
equipment.  Tables  4-17  and  4-18  list  the  maximum  AFs  and  their  corresponding 
frequencies  for  the  dry  and  wet  conditions,  respectively. 


Table  4-15 

SA  at  Frequency  of  50  Hz  for  Erection/Service  Bay  Block-Frame- 
Shear  Wall  Models,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g's 

SA  /  PGA 
g's 

SA50Hz 

g's 

SA  /  PGA 
g's 

14S9ErD110FW5 

110 

0.63 

3.32 

0.41 

2.56 

20S9ErD75FW5 

75 

0.42 

2.21 

0.31 

1.94 

12S9ErD45FW5 

45 

0.36 

1.89 

0.19 

1.19 

8S9ErD20FW5 

20 

0.2 

1.05 

0.17 

1.06 

Table  4-16 

SA  at  Frequency  of  50  Hz  for  Erection/Service  Bay  Block-Frame- 
Shear  Wall  Models,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g's 

SA  /  PGA 
g's 

SA50HZ 

g's 

SA  /  PGA 
g's 

14S9ErW110FW5 

110 

0.6 

3.16 

0.44 

2.75 

14S9ErW75FW5 

75 

0.47 

2.47 

0.32 

2.00 

12S9ErW45FW5 

45 

0.3 

1.58 

0.21 

1.31 

8S9ErW20FW5 

20 

0.25 

1.32 

0.17 

1.06 

Table  4-17 

Maximum  AF  and  Its  Corresponding  SAfor  Erection/Service  Bay 
Block-Frame-Shear  Wall  Models,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g's 

(AF)sa 
at  f 

f 

Hz 

SA 

g’s 

(AF)sa 
at  f 

14S9ErD110FW5 

110 

11 

2.41 

5.40 

10.37 

2.34 

7.13 

20S9ErD75FW5 

75 

15 

2.16 

6.36 

16 

1.51 

5.46 

12S9ErD45FW5 

45 

25 

1.11 

4.69 

27.6 

3.67 

2.02 

8S9ErD20FW5 

20 

- 

- 

- 

- 

- 

- 
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Table  4-18 

Maximum  AF  and  its  Corresponding  SAfor  Erection/Service  Bay 
Block-Frame-Shear  Wall  Models,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

f 

Hz 

SA 

g's 

(AF)sa 
at  f 

f 

Hz 

SA 

g's 

(AF)sa 
at  f 

14S9ErW110FW5 

110 

8.5 

3.51 

8.18 

8.38 

1.97 

5.78 

20S9ErW75FW5 

75 

13 

2.72 

6.64 

13 

1.9 

6.48 

12S9ErW45FW5 

45 

20 

1.30 

4.50 

20.89 

6.07 

2.77 

8S9ErW20FW5 

20 

33 

0.36 

1.68 

33 

0.19 

1.03 

The  PNW  maximum  AFs  are  less  than  those  for  the  CEUS  ground  motion  for 
the  wet  cases.  The  AF  curves  for  each  in  dry  and  wet  condition  models  are  pre¬ 
sented  in  Figures  4-18  through  4-25.  The  figures  were  obtained  for  both  ground 
motions  and  assuming  5  percent  damping  for  the  substructure.  The  AF  curves 
were  computed  by  using  response  spectrum  curves  with  5  percent  damping.  The 
results  for  2  percent  damping  are  presented  in  Appendix  F  for  use  with  the  equip¬ 
ment.  Table  4-19  lists  the  figure  numbers  for  the  AF  curves  obtained  for  the 
block- frame-shear  wall  models  of  the  erection/service  bay. 


Table  4-19 

Figure  Numbers  of  Erection/Service  Bay  Block-Frame-Shear  Wall 
Model  Results  for  5  Percent  Damping 

Height 

Ground  Motion  ! 

Group 

ft 

CEUS 

PNW 

1 

110 

Figure  4-18 

Figure  4-22 

2 

75 

Figure  4-19 

Figure  4-23 

3 

45 

Figure  4-20 

Figure  4-24 

4 

20 

Figure  4-21 

Figure  4-25 

4.7.2  Response  amplification  sensitivity  to  shear  wall  thickness  for 
the  erection/service  bay  block-frame-shear  wall  finite  element 
models 

The  erection/service  bay  block- frame-shear  wall  finite  element  model  dis¬ 
cussed  in  the  previous  section  includes  a  shear  wall  that  provides  lateral  resist¬ 
ance  to  the  ground  motion.  The  horizontal  ground  motion  was  applied  in  the 
upstream/downstream  direction  in  this  study.  The  shear  wall  thickness  was 
selected  by  inspection  of  the  construction  drawings  for  different  Corps  power¬ 
houses.  An  average  8-ft  wall  thickness  was  judged  to  represent  the  thickness  of 
the  shear  wall(s).  A  parametric  study  was  performed  in  order  to  investigate  the 
influence  of  the  shear  wall  thickness  on  the  response  amplification  of  the 
erection/service  bay  block- frame-shear  wall  model.  The  1 10-ft-high  substructure 
model  in  wet  condition  was  used  in  this  study.  The  thicknesses  of  the  shear  walls 
used  were  6,  8,  and  10  ft.  For  each  finite  element  model,  the  response  spectrum 
and  AF  at  top  of  the  substructure  were  computed. 


86 


Chapter  4  Time-History  Finite  Element  Results  for  the  Erection/Service  Bay  Monolith 


Amplification  Factor 


Frequency  (Hz) 
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c.  CEUS  ground  motion  response  spectrum 


Figure  4-18.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame 
shear  wall  model,  110-ft  height  (Dry  Model  ID:  14S9ErD110FW5, 
Wet  Model  ID:  14S9ErW1 10FW5),  for  CEUS  ground  motion  with 
5  percent  damping 
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Figure  4-19.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  75-ft  height  (Dry  Model  ID:  20S9ErD75FW5,  Wet 
Model  ID:  20S9ErW75FW5),  for  CEUS  ground  motion  with 
5  percent  damping 


Chapter  4  Time-History  Finite  Element  Results  for  the  Erection/Service  Bay  Monolith 


(6)  VS  (6)  VS 


Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 

1.6  -|  I  T  I  T  I 


1.4 

1.2 


1 


Frequency (Hz) 

.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


0.6  -r 

0.5  - 

0.4  - 

Ajt 

\ 

0.3  - 

/ 

V 

— 

0.2  - 

0.1  - 

1 

0  -1 

0  5  10  15  20  25  30  35  40  45  50 

Frequency (Hz) 
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Figure  4-20.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  45-ft  height  (Dry  Model  ID:  12S9ErD45FW5,  Wet 
Model  ID:  12S9ErW45FW5),  for  CEUS  ground  motion  with  5  percent 
damping 
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Figure  4-21.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  20-ft  height  (Dry  Model  ID:  8S9ErD20FW5,  Wet 
Model  ID:  8S9ErW20FW5),  for  CEUS  ground  motion  with  5  percent 
damping 
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Figure  4-22.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame 
shear  wall  model,  110-ft  height  (Dry  Model  ID:  14S9ErD110FW5, 
Wet  Model  ID:  14S9ErW1 10FW5),  for  PNW  ground  motion  with 
5  percent  damping 
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Figure  4-23.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  75-ft  height  (Dry  Model  ID:  20S9ErD75FW5,  Wet 
Model  ID:  20S9ErW75FW5),  for  PNW  ground  motion  with  5  percent 
damping 
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Figure  4-24.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame 
shear  wall  model,  45-ft  height  (Dry  Model  ID:  12S9ErD45FW5,  Wet 
Model  ID:  12S9ErW45FW5),  for  PNW  ground  motion  with  5  percent 
damping 
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Figure  4-25.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  20-ft  height,  (Dry  Model  ID:  8S9ErD20FW5,  Wet 
Model  ID:  8S9ErW20FW5),  for  PNW  ground  motion  with  5  percent 
damping 
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Table  4-20  lists  the  first  mode  period  and  frequency  for  the  three  1 10-ft- 
height  erection/service  bay  block- frame-shear  wall  models.  The  fundamental 
mode  frequency  increases  by  0.74  Hz  for  the  finite  element  models  with  an 
increase  in  shear  wall  thickness  from  6  to  10  ft.  Table  4-21  lists  the  values  of  the 
SA  at  50  Hz  (SA50hz)  and  its  ratio  with  the  PGA,  which  corresponds  to  the  AF  at 
a  frequency  of  50  Hz  at  top  of  the  substructure.  Table  4-22  lists  the  maximum  AF 
and  its  corresponding  frequency.  For  frequencies  near  the  first  natural  mode  of 
vibration  of  the  finite  element  models,  the  1 0-ft-thick  shear  wall  results  in  higher 
amplification  than  do  the  other  two  shear  wall  thicknesses. 


Table  4-20 

First  Mode  of  Vibration  and  PMR  for  the  110-ft-Height  Erection/ 
Service  Bay  Block-Frame-Shear  Wall  Models 

Model  ID 

Height 

ft 

First  Mode 

Final  Mode 

T i 

sec 

u 

Hz 

Total  No. 

/ 

PMR 

% 

T, 

sec 

fi 

Hz 

14S9ErD110FW5-6  ft 

110 

0.1252 

7.98 

18 

93.14 

0.014 

70.76 

14S9ErD110FW5-8  ft 

110 

0.1189 

8.41 

18 

93.3 

0.014 

71.11 

14S9ErD110FW5-10  ft 

110 

0.1147 

8.72 

18 

93.4 

0.014 

71.11 

Table  4-21 

SA  at  50  Hz  of  the  Three  Finite  Element  Models  Using  the  CEUS 
and  PNW  Ground  Motions 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g's 

SA50Hz/PGA 

g's 

SA50HZ 

g's 

sa50Hz/pga 

g's 

14S9ErD110FW5-6  ft 

110 

0.64 

3.37 

0.41 

2.56 

14S9ErD110FW5-8  ft 

110 

0.6 

3.16 

0.45 

2.81 

14S9ErD110FW5-10  ft 

110 

0.67 

3.53 

0.46 

2.88 

Table  4-22 

Maximum  AF  and  Its  Corresponding  Frequency  for  the  Three 
110-ft-Height  Erection/Service  Bay  Models _ 


Model  ID 

Height 

ft 

Maximum  AF 

CEUS 

PNW  | 

f,  Hz 

SA,  g's 

(AF)sa 
at  f 

f,  Hz 

SA,  g's 

(AF)sa 
at  f 

14S9ErD110FW5-6  ft 

110 

8 

4.09 

8.72 

8 

1.76 

5.62 

14S9ErD110FW5-8  ft 

110 

8.5 

3.59 

8.35 

8.5 

2.03 

5.82 

14S9ErD110FW5-10  ft 

110 

9 

4.69 

9.50 

9 

2.44 

6.83 

Because  the  ground  motion  response  spectra  used  for  the  analysis  are  not 
smooth-shaped,  broad-band  spectra,  a  small  change  in  frequency  may  result  in  an 
abrupt  change  in  results.  Figure  4-26  shows  the  SA  for  the  CEUS  and  PNW 
ground  motion  for  the  three  finite  element  models,  and  Figure  4-27  shows  the 
corresponding  AF  curves. 
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6)  VS  (6)  VS 


Figure  4-26.  Response  spectrum  curves  at  top-of-substructure  for  three  different 
shear  wall  thicknesses  of  the  erection  bay,  1 1 0-ft  height,  block- 
frame-shear  wall  model,  using  5  percent  damping 
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Figure  4-27.  AF  curves  at  top-of-substructure  for  three  different  shear  wall 

thicknesses  of  the  erection  bay,  110-ft  height,  block-frame-shear 
wall  model,  using  5  percent  damping 
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4.7.3  Response  amplification  sensitivity  to  location  at  top  of  the 
substructure  for  the  erection/service  bay  block-frame-shear  wall 
model 

Computed  results  reported  in  previous  sections  for  the  erection/service  bay 
block- frame-shear  wall  models  were  based  on  the  calculation  of  the  spectral 
acceleration  at  the  midnode  along  the  top  of  the  substructure.  In  order  to  investi¬ 
gate  the  variation  in  computed  results  along  the  top  of  the  substructure,  the  spec¬ 
tral  acceleration  and  AFs  were  obtained  at  two  other  locations  along  the  top-of- 
substructure.  For  this  parametric  study,  the  1 1 0-ft-height  block- frame-shear  wall 
model  in  the  wet  condition  was  used. 

Table  4-23  lists  the  values  of  the  SA  at  50  Hz  and  its  ratio  with  the  PGA, 
which  corresponds  to  the  AF  at  a  frequency  of  50  Hz.  Table  4-24  lists  the 
maximum  AF  and  its  corresponding  frequency  at  the  three  nodes  identified  in 
Figure  4-28.  Figure  4-29  shows  the  spectral  acceleration  for  the  three  locations 
along  the  top  of  the  substructure  using  the  CEUS  and  the  PNW  ground  motions. 
Figure  4-30  presents  the  AF  curves  for  the  three  different  nodes  along  the  top-of- 
substructure. 

The  results  show  that  the  response  of  the  substructure  is  higher  for  the 
upstream  node  location.  For  the  midnode  and  the  downstream  location  there  is  no 
significant  difference  in  response.  The  higher  response  at  the  upstream  node  loca¬ 
tion  is  due  to  the  added  water  mass  applied  to  the  upstream  face  of  the  substruc¬ 
ture,  as  explained  in  Appendix  D. 


Table  4-23 

SA  at  50  Hz  at  Three  Different  Locations  along  the  Top  of  the  110-ft-Height  Block- 
Frame-Shear  Wall  Substructure  for  the  Wet  Condition 

CEUS 

PNW 

Height 

SA50HZ 

SA50HZ 

Model  ID 

ft 

g's 

iSHHH 

g's 

iSttHHS 

14S9ErD110FW5-Upstream  node 

110 

0.76 

4 

0.51 

3.18 

14S9ErD110FW5-Midnode 

110 

0.6 

3.16 

0.45 

2.81 

14S9ErD110FW5-Downstream  node 

110 

0.69 

3.63 

0.46 

2.88 

Table  4-24 

Maximum  AF  and  Its  Corresponding  Frequency  at  Three  Different  Locations  along  the 
Top  of  the  110-ft-Height  Block-Frame-Shear  Wall  Substructure  for  the  Wet  Condition 


Model  ID 

Height 

ft 

Maximum  AF 

CEUS 

PNW  I 

f,  Hz 

SA,  g's 

(AF)sa 
at  f 

f,  Hz 

SA,  g's 

(AF)sa 
at  f 

1 4S9ErD1 1 0FW5-Upstream  node 

110 

8.5 

4.50 

9.96 

9 

2.55 

7.03 

1 4S9ErD1 1 0FW5-Midnode 

110 

8.5 

3.59 

8.35 

8.5 

2.03 

5.82 

14S9ErD110FW5-Downstream  node 

110 

8.5 

3.61 

8.39 

8.5 

2.04 

5.85 
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52  X-Y  Plane  @  Z=0 


Figure  4-28.  Location  of  nodes  in  the  1 1 0-ft-height  erection  bay  block-frame-shear  wall  model  used  in  the 
AF  calculations 
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Figure  4-29.  Response  spectrum  curves  at  three  different  locations  along  the  top- 
of-substructure  of  the  erection  bay,  1 1 0-ft-height  block-frame-shear 
wall  model,  in  wet  condition,  using  5  percent  damping 
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Figure  4-30.  AF  curves  at  three  different  locations  along  the  top-of-substructure 
of  the  erection  bay,  1 10-ft-height  block-frame-shear  wall  model,  in 
wet  condition,  using  5  percent  damping 
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5  Amplification  Factors  for 
Response  Spectra/Modal 
Analysis 


5.1.  Introduction 

The  results  from  time -history,  modal  analyses  of  powerhouse  substructures 
and  corresponding  AFs  are  presented  in  Chapter  3  for  the  generator  bay  substruc¬ 
ture  and  in  Chapter  4  for  the  erection/service  bay  substructure.  In  this  chapter,  the 
AF  curves  for  the  CEUS  and  for  the  PNW  ground  motions  are  plotted  in  the 
same  graph  but  distinguished  by  dry  and  wet  conditions. 

Figure  5-1  shows  the  response  spectra  at  rock  for  the  two  ground  motions 
considered  in  this  study.  The  two  SAs  were  normalized  by  their  respective  peak 
ground  accelerations  at  rock.  The  PNW  ground  motion  possesses  the  largest  SA 
value  of  2.8  g  at  a  frequency  of  4  Hz  compared  with  a  maximum  SA  value  of 
2.6  g  at  9  Hz  for  the  CEUS  ground  motions. 

In  Section  5.2,  the  computed  AF  curves  for  the  generator  bay  substructure 
are  discussed.  The  AF  curves  for  the  erection/service  bay  block  models  are  pre¬ 
sented  in  Section  5.3  and  for  the  block-frame-shear  wall  models  in  Section  5.4. 
An  example  comparison  between  (a)  a  complete  model  time-history  analysis, 

(b)  a  complete  model  response  spectra  analysis,  and  (c)  a  decoupled  response 
spectra  analysis  is  presented  in  Section  5.5  for  an  erection/service  bay  monolith. 
Finally,  in  Section  5.6,  a  demonstration  example  of  the  construction  of  a  top-of- 
substructure  PSA  versus  frequency  curve  is  given  for  an  idealized  erection/ 
service  bay  of  the  Cordell  Hull  Powerhouse  superstructure  using  a  National 
Earthquake  Hazards  Reduction  Program  (NEHRP)  (FEMA  2001)  design 
response  spectra  to  represent  the  in-rock  seismic  event. 

The  differences  between  SA  and  PSA  are  explained  in  Appendix  E.  In  that 
appendix  it  is  proved  that  for  a  low  percent  damping  and  short-period  structure, 
as  for  the  case  of  powerhouse  substructures,  the  differences  between  SA  and  PSA 
are  insignificant.  The  exact  solution  of  a  decoupled  analysis  of  an  isolated  super¬ 
structure  is  to  subject  it  with  a  total  acceleration  time-history.  Because  in  most  of 
the  cases,  only  the  maximum  response  is  of  concern  and  not  the  response  of  the 
structure  in  time,  the  ground  motion  is  approximated  by  a  spectral  acceleration 
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Figure  5-1 .  In-rock  CEUS  and  PNW  normalized  response  spectra  versus 
frequency  (5  percent  damping) 

and  a  modal/response  spectra  analysis  is  conducted.  For  simplicity  in  the  calcula¬ 
tions,  and  because  of  a  low  damping  and  short-period  structure,  a  PSA  is  com¬ 
monly  used  as  input  at  the  base  of  the  finite  element  models.  The  PSA  is  then 
transferred  from  the  rock  to  the  top-of-substructure  by  the  use  of  the  AF  curves 
developed  in  this  study.  The  AF  curves  were  calculated  by  using  absolute  accel¬ 
erations  that  take  into  account  the  elastic  and  damping  effects  of  the  substructure, 
and  are  then  approximately  transferred  to  the  superstructure. 


5.2.  Computed  AF  Relationships  for  Generator 
Bay  Composite  Models 

The  AF  curves  in  Figures  5-2  to  5-5  for  the  generator  bay  substructure  were 
obtained  for  the  CEUS  and  PNW  ground  motions  for  dry  and  wet  conditions.  A 
value  of  5  percent  damping  is  assumed  appropriate  for  the  substructure  in  all 
dynamic  analyses.  The  response  spectra  at  rock  and  at  top-of-substructure  used 
for  the  development  of  the  AF  curves  were  developed  using  5  percent  damping. 
This  value  of  damping  is  typically  recommended  for  dynamic  response  spectra, 
modal  analysis  of  a  superstructure  to  earthquake  shaking.  The  AF  curves  for 
analyzing  equipment  resting  on  top-of-substructure,  which  were  developed  using 
response  spectra  with  2  percent  damping,  are  presented  in  Appendix  F. 
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Figure  5-3.  AF  versus  frequency  curves,  generator  bay,  composite  model,  100-ft 
height  with  5  percent  damping 
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Figure  5-4.  AF  versus  frequency  curves,  generator  bay,  composite  model,  70-ft 
height  with  5  percent  damping 
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Figure  5-5.  AF  versus  frequency  curves,  generator  bay,  composite  model,  40-ft 
height  with  5  percent  damping 
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5.3.  Computed  AF  Relationships  for  Erection/ 
Service  Bay  Block  Models 

The  AF  curves  presented  in  Figures  5-6  to  5-9  are  for  the  erection/service 
bay  block  models.  They  were  obtained  using  5  percent  damping  and  response 
spectra  for  the  acceleration  time-histories  at  rock  and  at  the  top-of-substructure. 
For  dynamic  response  spectra,  modal  analysis  of  a  superstructure  to  earthquake 
shaking,  the  5  percent  damping  value  is  recommended.  The  AF  curves  for 
analyzing  equipment,  which  were  developed  using  response  spectra  with  2  per¬ 
cent  damping,  are  presented  in  Appendix  F.  For  the  dynamic  analysis  of  the  sub¬ 
structure,  5  percent  damping  was  assumed  in  the  finite  element  block  models. 


5.4.  Computed  AF  Relationships  for  Erection/ 
Service  Bay  Block-Frame-Shear  Wall  Model 

A  second  type  of  finite  element  model  that  was  used  to  idealize  the  erection/ 
service  bay  substructure  is  referred  to  as  a  block- frame-shear  wall  model.  It 
consists  of  a  U-ffame  shaped  mass  concrete  structure  with  structural  frames  and 
shear  walls  inside  the  U  of  the  model.  The  geometry  and  details  of  the  models 
and  discussion  of  the  AF  were  presented  in  Chapter  4.  This  section  summarizes 
the  AF  curves  for  the  CEUS  and  PNW  ground  motions. 

Figures  5-10  to  5-13  show  the  AF  curves  for  each  of  the  two  ground  motions. 
Results  were  divided  for  each  block- frame-shear  wall  model  group  into  dry  and 
wet  conditions  and  were  developed  using  5  percent  damping. 


5.5  Comparison  Example:  Time-History/Modal 
Response,  Response  Spectra  Analysis,  and 
Decoupled  Response  Spectra  Analysis  of  an 
Erection/Service  Bay 

This  section  discusses  a  comparison  of  the  results  of  a  complete  dynamic 
response  spectra  analysis  and  a  decoupled  response  spectra  analysis.  The  finite 
element  models  considered  are  of  an  erection/service  bay  substructure  and  super¬ 
structure.  The  superstructure  dimensions  were  an  idealization  of  the  Cordell  Hull 
Powerhouse  that  was  in  Group  2,  the  75-ft- height  block  model,  assuming  a  wet 
condition.  The  ground  motion  considered  was  the  representative  one  of  the 
CEUS.  A  value  of  5  percent  damping  was  used  for  the  response  spectra  at  rock 
and  at  top-of-substructure  calculations. 

The  structural  analysis  software  SAP2000  was  used  to  create  and  analyze  the 
finite  element  models.  Two  models  were  constructed  for  this  example.  The  first 
model  is  identified  as  a  complete  model.  It  consists  of  the  75-ft-height  block 
model  for  the  substructure  with  a  frame  structure  constructed  on  top  of  a  solid 
“block”  of  mass  concrete.  The  Cordell  Hull  superstructure  was  used  as  a  guide 
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(Amplification  Factor)SA  (Amplification  Factor)SA 
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Figure  5-7.  AF  versus  frequency  curves,  erection/service  bay,  block  model,  75-ft 
height  with  5  percent  damping 
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Figure  5-8.  AF  versus  frequency  curves,  erection/service  bay,  block  model,  45-ft 
height  with  5  percent  damping 
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Figure  5-9.  AF  versus  frequency  curves,  erection/service  bay,  block  model,  20-ft 
height  with  5  percent  damping 
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Figure  5-10.  AF  versus  frequency  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  1 1 0-ft  height  with  5  percent  damping 
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- CEUS 

. PNW 


Frequency  (Hz) 

b.  Wet  condition 


Figure  5-11.  AF  versus  frequency  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  75-ft  height  with  5  percent  damping 
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Figure  5-12.  AF  versus  frequency  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  45-ft  height  with  5  percent  damping 
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(Amplification  FactorJsA 


Figure  5-13.  AF  versus  frequency  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  20-ft  height  with  5  percent  damping 
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for  the  properties  of  the  frame  superstructure.  The  second  model  is  identified  as 
the  decoupled  frame  model.  It  consists  of  only  the  frame  that  represents  the 
superstructure.  The  structure  material  properties  assumed  for  this  example  was 
mass  concrete  and  reinforced  concrete  with  the  properties  listed  in  Table  5-1  for 
the  mass  of  concrete  and  Table  5-2  for  the  frames.  Note  that  the  type  of  finite 
elements  and  their  material  properties  were  previously  presented  in  Table  4-2  for 
the  block  model.  Table  5-3  lists  the  idealized  geometry  of  the  elements  used  in 


SAP2000. 


Table  5-1 

Material  Properties  Used  in  the  Comparative  Example  for  the  Plane 

Finite  Elements 

|  Mass  Concrete  Properties  | 

Compressive  strength  f  'c 

3,750 

psi 

Modulus  of  elasticity  Ec 

3.7E+6 

psi 

Mass  of  concrete  mc 

4.6584 

Ib*sec7ft 

Unit  weight  yc 

150 

pcf 

Poisson’s  ratio  u 

0.2  i 

Table  5-2 

Material  Properties  Used  in  the  Frame  Elements 

Mass  Concrete  Properties 

Compressive  strength  f  'c 

4,500  psi 

Modulus  of  elasticity  Ec 

4.1  E+6  psi 

Mass  of  concrete  mc 

4.6584  Ib*sec7ft 

Unit  weight  yc 

150  pcf 

Poisson’s  ratio  u 

0.2 

Table  5-3 

Geometry  of  Comparative  Example  SAP2000  Model 

Element 

Dimensions 

ft 

9-node  plane  stress  elements,  thickness 

43 

Substructure  vertical  frame  elements 

1.5  by  43 

Substructure  horizontal  frame  elements 

1  by  43 

Superstructure  horizontal  frame  element 

0.10  by  43 

Superstructure  top  vertical  frame  element 

2.5  by  43 

Superstructure  bottom  vertical  frame  element 

5  by  43 

The  “complete”  finite  element  model  is  shown  in  Figure  5-14.  The  response 
spectrum  and  time-history  at  rock  used  in  the  two  complete  analyses  are  also 
shown  in  this  figure.  The  earthquake  used  in  the  analysis  is  the  CEUS  ground 
motion  for  the  development  of  the  AF  curves  in  this  study.  Figure  5- 15a  shows 
the  decoupled  frame  model.  It  consists  of  a  single  bay  frame  structure  modeled 
with  9  nodes,  each  containing  lumped  masses  in  each  vertical  element.  The 
decoupled  frame  model  was  subjected  to  the  amplified  response  spectrum  shown 
in  Figure  5-1 5b  that  was  constructed  using  the  AFs  obtained  from  the  75-ft- 
height  block  model  assuming  wet  conditions.  The  response  spectrum  in 
Figure  5- 19b  was  developed  by  multiplying  the  in-rock  CEUS  SA  value  at  each 
frequency  by  the  AF  for  the  same  frequency. 
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a.  SAP2000  finite  element  model-complete  analysis 


b.  Rock  response  spectra  for  complete  analysis  (CEUS  ground  motion  -  5  percent  damping) 


c.  Rock  time-history  acceleration  for  complete  analysis  (CEUS  ground  motion) 


Figure  5-14.  SAP2000  complete  response  spectra  and  time-history  analysis  of 
Cordell  Hull  erection/service  bay  cross  section  for  comparison 
example 
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lane  @ Z=0  |  X38.08  Y54  84  ZO  00  [km!  7| 


a.  SAP2000  frame  superstructure  for  decoupled  analysis 


b.  CEUS  amplified  pseudo-response  spectrum  (top-of-substructure)  using  the  75-ft-height  block  model  AFs,  wet  condition 


Figure  5-15.  SAP2000  decoupled  analysis  of  Cordell  Hull  Powerhouse  for  comparative  example 
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Table  5-4  lists  the  maximum  displacements  computed  at  the  top  of  the  frame 
superstructure.  There  is  not  a  significant  difference  in  maximum  displacement 
between  the  results  of  the  three  analyses.  Tables  5-5  and  5-6  present  the  internal 
forces  calculated  at  different  points  in  the  frame  superstructure.  The  locations  of 
the  points  in  Tables  5-5  and  5-6  are  shown  in  Figure  5-16.  In  these  tables  also  are 
listed  the  percent  of  difference  between  the  response  spectra  analysis  of  the  com¬ 
plete  model  and  the  decoupled  frame  model  with  the  amplified  response  spectra. 
The  maximum  difference  in  shear  forces  was  17.6  percent.  It  was  found  in  the 
left  side  of  the  crane  at  support  point  2.  A  maximum  difference  of  17.1  percent  in 
moment  forces  was  calculated  at  the  right  side  in  the  crane  (point  4). 


Table  5-4 

Maximum  Displacement  at  Top  of  Frame  Superstructure 

Model  ID 

Description 

CEUS 

ft 

PNW 

ft 

CordellHullW-Complete_TH 

Complete  Model 

Time-History  Analysis 

0.04 

0.05 

CordellHullW-Complete_RS 

Complete  Model 

Response  Spectra  Analysis 

0.044 

0.053 

CordellHullW-Frame_RS 

Decoupled  Frame  Model 

Block  Response  Spectrum 

0.0419 

0.0511 

Table  5-5 

Shear  Forces  Within  Frame  Superstructure  for  CEUS  Earthquake 

Point 

Model 

CompleteRS 

kip 

Decoupled  Frame  RS 

CompleteTH 

kip 

kip 

%  Difference 

1  (LHS-Base) 

439 

385.5 

12.2 

-445 

2  (LHS-Crane) 

334.4 

275.65 

17.6 

-171.5 

3  (RHS-Base) 

501.3 

456.2 

9.0 

-576.8 

4  (RHS-Crane) 

313 

290.4 

7.2 

-167.8 

Table  5-6 

Moment  Forces  Within  Frame  Superstructure  for  CEUS  Earthquake 

Point 

Model 

CompleteRS 

kip-ft 

Decoupled  Frame  RS 

CompleteTH 

kip-ft 

kip-ft 

%  Difference 

1  (LHS-Base) 

10,042 

9,351.9 

6.9 

9,993 

2  (LHS-Crane) 

3,614.03 

3,029.9 

16.2 

3,983 

3  (RHS-Base) 

10,689 

10,155.4 

5.0 

8,733 

4  (RHS-Crane) 

3,770.64 

3,124.09 

17.1 

3,423 
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Figure  5-16.  Frame  lumped  mass  model 

5.6.  Demonstration  Example:  Decoupled  Analysis 
of  an  Erection/Service  Bay  Using  NEHRP  Maps 

In  this  section,  the  procedure  for  a  decoupled  dynamic  analysis  is  outlined.  A 
superstructure,  represented  by  a  frame  structure,  is  analyzed  by  a  response 
spectra/modal  analysis  procedure  using  a  design  response  spectrum  developed  by 
following  the  NEHRP  procedure  (FEMA  2001).  The  in-rock,  design  response 
spectrum  was  amplified  by  the  AF  developed  in  this  study  to  obtain  the  approxi¬ 
mate  top-of-substructure  response  spectrum  for  the  frame  superstructure.  The 
computed  maximum  displacements  at  top  of  the  frame,  shear  forces,  and  moment 
forces  are  reported. 

A  frame  structure  located  at  the  top  of  an  erection/service  bay  substructure 
was  analyzed  for  earthquake  loadings.  The  powerhouse  is  located  at  latitude  of 
36.29  and  longitude  of -85. 9416  in  the  central-eastern  part  of  the  United  States. 

It  was  assumed  that  the  frame  structure  was  constructed  with  reinforced  concrete. 
Its  properties  are  listed  in  Table  5-2.  From  construction  drawings,  an  estimation 
of  the  dimensions  of  the  structural  elements  was  made,  resulting  in  the 
simplified/idealized  dimensions  listed  in  Table  5-7.  Figure  5-17  shows  a  sche¬ 
matic  of  the  frame  that  was  analyzed.  The  superstructure  horizontal  element  is 
identified  as  SH,  the  superstructure  top  vertical  element  is  identified  as  STV,  and 
the  superstructure  bottom  vertical  element  as  SBV. 

The  in-rock,  design  response  spectrum  used  in  this  analysis  corresponds  to  a 
2,475-year  return  period  and  was  developed  for  the  powerhouse  site  according  to 
the  information  obtained  at  http://earthqnake.usgs.gov/research/hazmaps/design/. 
From  this  map  the  peak  ground  acceleration  at  rock  (PGArock)  is  0. 1  lg,  the  PSA 
at  a  period  of  0.2  sec  (PSA0.2)  is  0.27g,  and  the  PSA  at  a  period  of  1.0  sec  (PSAi) 
is  0.1 1 18g.  Figure  5-18  shows  the  resulting  in-rock  design  response  spectrum. 
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Table  5-7 

Geometry  of  the  Frame  Structural  Elements  for  Demonstration 
Example 

Element 

Description 

Dimensions 

ft 

SH 

Superstructure  horizontal  frame  element 

0.10  by  43 

STV 

Superstructure  top  vertical  frame  element 

2.5  by  43 

SBV 

Superstructure  bottom  vertical  frame  element 

5  by  43 

Figure  5.17.  Frame  superstructure  for  demonstration  example 

The  design  response  spectrum  at  rock,  determined  by  NEHRP  procedure,  is 
transferred  to  the  top  of  the  substructure  for  the  decoupled  dynamic  analysis. 
Assuming  a  wet  condition,  the  AF  curve  for  the  75-ft  block  model  of  the 
erection/service  bay  substructure  (Group  2)  is  used.  The  PSA  at  rock  at  each 
frequency  is  multiplied  by  the  AF  at  that  same  frequency.  The  resulting  top-of- 
substructure  amplified  design  response  spectrum  is  shown  in  Figure  5-19. 

The  structural  analysis  software  SAP2000  was  then  used  for  creating  a 
simplified/idealized  finite  element  model  and  for  the  dynamic  analysis  of  the 
superstructure  only.  The  model  consists  of  frame  elements  with  three  degrees  at 
each  node  (two  translations  and  one  rotation)  with  nine  lumped  masses  specified 
along  its  vertical  elements.  A  response  spectrum/modal  analysis  was  then  per¬ 
formed  to  calculate  the  maximum  displacement  at  the  top  of  the  frame  super¬ 
structure  and  the  shear  and  moment  forces  at  the  base  and  at  the  crane  levels. 
Figure  5-20  shows  the  SAP2000  model  used  in  this  example. 
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Figure  5-18.  Design  rock  response  spectrum  constructed  using  the  NEHRP  procedure  with  PGA  anchor 
from  USGS  map  for  demonstration  example  site 
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Figure  5-19.  Amplified  USGS/NEHRP  response  spectrum  at  top-of-substructure 
constructed  using  the  75-ft  block  model  AF  curve  for  the  wet 
condition 


lane  @  Z-0  X38.08  Y54.84  ZO.OO  “IIS  S 

Figure  5-20.  SAP2000  frame  superstructure  lumped  mass  model 
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The  analysis  results  are  summarized  in  Table  5-8.  A  maximum  displacement 
of  0.0366  ft  at  top  of  the  frame  was  calculated  using  the  NEHRP  procedure 
with  the  75-ft  block  AF.  The  maximum  shear  force  calculated  for  the  right-side 
column  at  the  base  point  3  was  3 16.37  kips  with  a  maximum  moment  of 
7,926  kip-ft. 


Table  5-8 

Shear  Forces  and  Moments  within  Frame  from  an  NEHRP-Based, 
Amplified  Design  Response  Spectrum 

Shear 

Moments 

NEHRP 

NEHRP 

Block  AF 

Block  AF 

Point 

kip 

kip-ft 

1  (LHS-Base) 

282.48 

7,677 

2  (LHS-Crane) 

214.47 

2,455.99 

3  (RHS-Base) 

316.37 

7,926.14 

4  (RHS-Crane) 

219.42 

2,938.64 
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6  Summary  and  Conclusions 


This  report  summarizes  the  results  of  a  series  time-history  analysis  of 
powerhouse  substructures  conducted  in  order  to  develop  amplification  factor 
(AF)  curves  (i.e.,  AF  versus  frequency)  at  the  top  of  the  substructure  to  facilitate 
a  decoupled  response  spectrum/modal  analysis  of  the  superstructure  (or 
equipment).  The  massive  size  of  the  substructure  allows  for  this  decoupled 
dynamic  analysis  of  the  superstructure  and/or  equipment.  In  this  simplified 
procedure  of  modal  analysis,  AF  curves  are  used  to  transfer  the  rock  response 
spectra  to  the  top  of  the  substructure  for  use  in  the  decoupled  analysis. 

The  construction  of  the  amplification  curves  used  in  the  decoupled  response 
spectra/modal  analysis  of  a  superstructure  or  equipment  consisted  of  three  steps: 

•  For  a  specified  in-rock  acceleration  time-history,  construct  a  pseudo¬ 
response  spectra. 

•  Construct  a  finite  element  model  of  the  substructure  only  and  using 
SAP2000  compute  the  top-of-substructure  acceleration  time-history 
response. 

•  Construct  the  AF  versus  frequency  curve  by  computing  the  ratio  of  the 
top-of-substructure  value  for  absolute  acceleration  divided  by  the  in-rock 
value  absolute  acceleration  for  each  frequency  up  to  50  Hz. 

Two  acceleration  time-histories  and  their  corresponding  response  spectra 
were  developed  representing  an  earthquake  event  in  the  central-eastern  part  of 
the  United  States  (CEUS)  and  an  event  in  the  Pacific  Northwest  of  the  United 
States  (PNW).  They  were  developed  synthetically  for  a  seismic  hazard  exposure 
period  of  1,000  years.  The  peak  ground  acceleration  at  rock  for  the  CEUS 
earthquake  was  0.1 9g  with  a  maximum  pseudo-spectral  acceleration  of  0.49g  at 
frequency  of  9  Hz  (5  percent  damping).  For  the  PNW  earthquake  the  maximum 
acceleration  was  0. 16g  with  maximum  pseudo-spectral  acceleration  of  0.46g  at  4 
Hz  (5  percent  damping). 

A  “box”  model  consisting  of  five  structural  members  idealized  the  generator 
bay  substructure  monoliths.  Two  shear  walls  along  the  upstream/downstream 
direction,  two  flexural  walls,  and  a  base  slab  constituted  the  idealized  geometry. 
Coips  generator  bay  substructures  were  tabulated,  evaluated,  and  generalized  into 
four  groups  according  to  their  heights.  Four  finite  element  models  were  created 
and  analyzed  using  the  capabilities  of  the  structural  analysis  software  SAP2000. 
The  analyses  were  performed  considering  the  normal  pool  and  tailrace  water 
elevations  (wet  condition)  and  for  the  dry  condition. 
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Two  models  were  created  to  idealize  the  erection/service  bay  substructure. 
One  model  consisted  of  a  block  finite  element  model,  and  the  second  a  block- 
frame-shear  wall  model.  A  block  model  represents  a  stiff  substructure,  and  the 
block- frame-shear  wall  model  represents  a  more  flexible  substructure  where  the 
inside  frames  influence  the  structural  system  response  to  earthquake  loadings. 

The  erection/service  bays  were  also  categorized  into  four  groups  according  to 
their  heights. 

The  AF  curves  at  top-of-substructure  for  each  model  were  calculated  using 
the  results  from  time-history  analysis  as  previously  described.  For  the  generator 
bay  models,  the  wet  conditions  result  in  greater  AFs  than  the  dry  cases  for  the 
125-ft-  and  40-ft-height  models  using  the  CEUS  ground  motion.  The  generator 
bay  of  group  4  (40-ft  height)  did  not  show  a  significant  amplification  for  the 
PNW  earthquake  assuming  dry  conditions.  The  maximum  AF  for  the  wet 
condition  was  found  for  group  1  (125-ft  height)  to  be  a  value  of  7.3  at  a 
frequency  of  8  Hz. 

The  AF  curves  were  also  developed  for  the  block  and  block- frame-shear  wall 
models  of  the  erection/service  bay.  The  maximum  AF  for  the  wet  condition  of 
the  block  models  is  6.95  at  a  frequency  of  9.63  Hz  for  the  PNW  ground  motion. 
For  the  block- frame-shear  wall  models,  the  maximum  AF  is  8. 1 8  at  a  frequency 
of  8.5  Hz  for  the  CEUS  ground  motion. 
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Appendix  A 
Summary  of  Generator 
and  Erection/Service  Bay 
Monolith  Geometry 


Fifty-four  Corps  powerhouse  plants  were  reviewed  in  this  study.  The  gen¬ 
erator  bay  and  erection  bay  monolith  cross  sections  were  evaluated  and  cate¬ 
gorized  into  one  of  four  groups  according  to  height  of  substructure.  Table  A-l 
lists  a  summary  of  the  heights  considered  and  the  height  assigned  to  each  of  the 
four  generator  and  four  erection/service  bay  groups.  Table  A-2  presents  a  list 
sorted  by  height  of  the  generator  bay  substructures  and  their  respective  widths. 
Table  A-3  presents  a  list  sorted  by  height  of  the  erection/service  bay  substructure 
and  their  dimensions. 


Table  A-1 

Summary  of  Heights  Considered  for  Each  Generator  and 
Erection/Service  Bay  Group _ 


Generator  Bay  i 

Height  H 

Number  of 

Maximum 

Minimum 

Assigned 

Hydropower 

Group  Number 

ft 

ft 

ft 

Plants1 

1 

135 

115 

125 

5 

2 

115 

85 

100 

19 

3 

85 

55 

70 

19 

4 

55 

- 

40 

11 

Erection/Service  Bay  j 

Height  H 

Number  of 

Maximum 

Minimum 

Assigned 

Hydropower 

Group  Number 

ft 

ft 

ft 

Plants1 

1 

135 

91 

110 

11 

2 

90 

60 

75 

13 

3 

55 

30 

45 

19 

4 

29 

- 

20 

7 

1  Not  all  Corps  hydropower  plants  were  included  in  this  evaluation. 
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Table  A-2 

Powerhouses  Sorted  by  Generator  Bay  Heights 

Group 

Number 

Generator  Bay  i 

Project 

Height 

H 

ft 

Upstream/Downstream 
Width  W 
ft 

Plan  Width 

B 

ft 

John  Day 

135 

145 

90 

Bonneville  II 

130 

130 

94 

3 

o 

Millers  Ferry 

120 

164 

80 

o 

Big  Bend 

120 

140 

86 

Gavins  Point 

118 

120 

80 

Bonneville  1 

113 

134 

62 

Dardanelle 

110 

123 

80 

Carters 

108 

115 

63 

Barkley 

108 

150 

86 

Cordell  Hull 

106 

184 

86 

The  Dalles 

106 

124 

86 

Ice  Harbor 

105 

138 

86 

Richard  B.  Russell 

100 

115 

71 

CM 

Walter  F.  George 

100 

118 

67 

3 

o 

Robert  F.  Henry 

100 

146 

73 

o 

Cheatham 

100 

169 

82 

Garrison 

96 

111 

75 

Harry  S.  Truman 

96 

156 

55 

Fort  Randall 

95 

102 

70 

West  Point 

90 

117 

77 

Jim  Woodruff 

90 

124 

67 

Chief  Joseph 

90 

102 

70 

Oahe 

90 

122 

76 

Sam  Rayburn 

90 

123 

55 

Allatoona 

75 

124 

54 

J.  Strom  Thurmond 

70 

93 

62 

Big  Cliff 

70 

72 

40 

Laurel 

70 

128 

78 

Center  Hill 

70 

137 

58 

Greers  Ferry 

70 

90 

56 

Table  Rock 

70 

131 

54 

Ozark 

69 

206 

65 

CO 

Lookout  Point 

68 

95 

60 

3 

o 

Webbers  Falls 

67 

206 

70 

o 

Fort  Peck 

65 

85 

60 

Wolf  Creek 

65 

137 

58 

J.  Percy  Priest 

65 

107 

80 

John  S.  Kerr 

63 

73 

70 

Dale  Hollow 

62 

115 

42 

Dexter 

60 

93 

90 

Sault  Ste.  Marie 

60 

74 

54 

Beaver 

60 

94 

61 

Blakely  Mountain 

60 

66 

56 

Buford 

54/64 

110 

62 

Fort  Gibson 

55 

80 

53 

Whitney 

55 

90 

46 

Bull  Shoals 

53 

79 

54 

■o- 

Broken  Bow 

52 

104 

53 

3 

o 

Norfork 

47 

58 

54 

<5 

Hills  Creek 

45 

55 

42.52 

Sault  Ste.  Marie  (old) 

40 

50 

49 

Foster 

35 

43 

40 

Philpott 

25 

40 

38 

Cougar 

25 

47 

34 

Appendix  A  Summary  of  Generator  and  Erection/Service  Bay  Monolith  Geometry 


Table  A-3 

Powerhouses  Sorted  by  Erection/Service  Bay  Heights 

Group 

Number 

Erection  or  Service  Bay  ] 

Project 

Height 

ft 

Upstream/Downstream 
Width,  ft 

Plan  Width 
ft 

Ozark 

- 

- 

65 

Sault  Ste.  Marie 

- 

74 

54 

Beaver 

- 

68 

65 

Fort  Gibson 

- 

- 

57 

Bonneville  II 

135 

150 

149 

3 

o 

Millers  Ferry 

120 

164 

80 

0 

Bonneville  1 

113 

134 

45 

Barkley 

108 

150 

86 

John  Day 

100 

145 

79/90 

The  Dalles 

100 

124 

86 

Garrison 

91 

111 

75 

West  Point 

90 

117 

90 

Jim  Woodruff 

90 

119 

42 

Sam  Rayburn 

90 

123 

45 

Robert  F.  Henry 

85 

146 

50 

Cordell  Hull 

80 

163 

43 

CM 

Big  Bend 

75 

140 

86 

3 

o 

Carters 

65 

115 

75 

0 

Fort  Peck 

65 

85 

60 

Dardanelle 

64 

125 

74 

Walter  F.  George 

63 

118 

55 

Cheatham 

62 

147 

60 

Greers  Ferry 

60 

90 

63 

John  S.  Kerr 

60 

73 

69 

Richard  B.  Russell 

55 

128 

100 

Chief  Joseph 

55 

148 

73 

Allatoona 

50 

110 

54 

Center  Hill 

50 

137 

64 

Oahe 

47 

97 

76 

Norfork 

47 

58 

54 

Table  Rock 

46 

114 

68 

Buford 

45 

110 

43 

CO 

Wolf  Creek 

45 

137 

66 

3 

o 

Dale  Hollow 

45 

115 

68 

0 

Bull  Shoals 

45 

79 

48 

Ice  Harbor 

44 

146 

83 

Whitney 

44 

90 

48 

Gavins  Point 

42 

120 

92 

Lookout  Point 

41 

102 

60 

Sault  Ste.  Marie  (old) 

40 

50 

49 

Blakely  Mountain 

36 

66 

56 

Fort  Randall 

34 

102 

70 

Harry  S.  Truman 

30 

156 

77 

Broken  Bow 

29 

94 

49 

Big  Cliff 

25(approx) 

72 

40 

Philpott 

25 

40 

20 

Hills  Creek 

20 

55 

39 

''T 

Foster 

0 

39 

40 

3 

o 

Cougar 

0 

47 

47 

0 

J.  Strom  Thurmond 

None 

Laurel 

None 

Webbers  Falls 

None 

J.  Percy  Priest 

None 

Dexter 

None 
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Appendix  B 

Equivalent  Flexural  Wall  in 
Generator  Bay  Idealized 
Geometry 


The  model  to  idealize  the  generator  bay  substructure  consisted  of  a  com¬ 
posite  model  in  two  dimensions  of  five  structural  features  that  form  a  box  as 
shown  in  Figure  B-la.  Figure  B-lb  shows  a  three-dimensional  view  of  the 
construction  of  the  two-dimensional  composite  model  used  in  the  finite  element 
analyses.  Shear  Wall  2  was  superimposed  on  shear  Wall  1  to  have  both  in  the 
same  plane.  The  open  space  in  the  center  of  the  box  (Figure  B-la)  represents  the 
location  of  the  turbine. 

The  flexural  walls  are  identified  by  Wall  4  facing  the  upstream  side  and 
Wall  5  facing  the  downstream  side.  The  real  generator  bay  substructure  upstream 
walls  have  the  geometry  shown  in  the  left  side  of  Figure  B-2.  To  simplify  the 
model,  the  two  walls  that  form  the  channel  to  the  water  intake  at  the  upstream 
face  were  converted  to  an  equivalent  rectangular  flexural  wall  (Wall  4).  Applying 
the  Parallel  Axis  Theorem,  the  inertia  of  the  real  geometry  was  calculated  at  the 
center  of  the  generator  bay,  and  an  equivalent  flexural  wall  area  was  calculated 
using  the  actual  (real)  wall  inertia.  Table  B-l  lists  a  summary  of  the  structural 
features  dimensions  for  the  generator  bay  substructure  model  for  Groups  1  and  2. 
In  this  appendix,  the  equivalent  inertia  calculations  are  presented  as  written  in 
Maple  6  code. 
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Figure  B-1 .  Three-dimensional  view  of  the  idealized  composite  structural  model  for  generator  bay 
substructure  idealization 
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Real  Generator  Bay  Monolith 
Plan  View 


Generator  Bay  Monolith  Idealization 
Plan  View 


Shear  Wall 


Real  Inertia 
Flexural  Wall 


Shear  Wall 


Flexural  Wall 


Shear  Wall 


4  |  a 

■eq»  ^eq 


Equivalent  Inertia  1  Flexural  Wall 

Flexural  Wall  shear  Wall 


Figure  B-2.  Plan  view  of  the  real  generator  bay  channel  walls  and  its  conversion  to  equivalent  inertia 
flexural  wall 
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Table  B-1 

Summary  of  Structural  Features  Dimensions  for  Generator  Bay  Substructure  Model 

Wall  4 

Wall  4 

Wall  5 

Wall  5 

Wall  1 

Wall  2 

Base  3 

Real 

Ideal 

Real 

Ideal 

Pool  AH1 

Tail  AH1 

Powerhouse 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

Group  1  (Height  =  125  ft,  Width  = 

135  ft,  Plan  Width  =  85  ft) 

John  Day 

Bonneville  II 

Millers  Ferry 

15 

15 

10 

(2)  7x53 

6 

(2)  6x57 

7 

2 

-35 

Big  Bend 

12 

15 

17 

(2)  7x72 

14 

(2)  5x57 

8 

37 

-32 

Gavins  Point 

13 

12 

10 

(2)  7x60 

10 

(2)  7x55 

8 

13 

-33.6 

Group  2  (Height  =  100  ft,  Width  = 

130  ft,  Plan  Width  =  75  ft) 

Bonneville  1 

25 

8 

16 

17 

-15 

Dardanelle 

12 

12 

10 

(2)  8  x  50 

8 

(2)  8x48 

7 

9 

-20 

Carters 

Barkley 

12 

12 

9 

(2)  7x52 

7 

(2) 6x50 

5 

Cordell  Hull 

15 

12 

10 

(2)  7x50 

9 

(2)  7x54 

8 

13 

-34 

The  Dalles 

Ice  Harbor 

80.75 

-14.25 

Richard  B.  Russell 

Walter  F.  George 

13 

12 

9 

45.5 

19 

Robert  F.  Henry 

10 

13 

16 

4 

Cheatham 

15 

12 

10 

(2)  7x50 

9 

(2)  6x57 

7 

-12 

-36 

Garrison 

10 

15 

15 

(2)  6x45 

5 

Harry  S.  Truman 

Fort  Randall 

15 

17 

10 

West  Point 

14 

12 

14 

(2)  5x58 

7 

(2)  5x58 

7 

38 

-25 

Jim  Woodruff 

12 

12 

6 

(1)6x36 

2 

(1)5x50 

1 

-3 

-36 

Chief  Joseph 

Oahe 

15 

15 

11 

Sam  Rayburn 

1  AH  as  measured  from  top  of  structure. 

Maple  6  Example  of  Calculations 

>  restart: 

Calculation  of  equivalent  moment  of  inertia  of  flexural  Wall  4  for 
generator  idealization. 

Inertia  for  actual  geometry 
Input 

bact :=8 : 
hact :=48 : 
dact :=55 : 

Actual  moment  of  inertia 

Iact :=2* ( ( (bact*hactA3) /12) + (bact*hact) *dactA2) ; 
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Inertia  for  flexural  wall  idealization 
Input 

W:=135 : 
bw : =85 : 

Ifw:=( (bw*hwA3) /3) - 

( (W*bw*hwA2 ) /2 ) + ( (bw*hw*WA2 ) / 4 ) =Iact 

Solve  for  hw: 

hw:=solve (Ifw,hw) ; 


hw : =  7.1000 


Equivalent  Flexural  Wall  in  Generator  Bay  Idealized  Geometry 


Appendix  C 

Investigation  of  Plane  Stress 
and  Plane  Strain  Modal  Time- 
History  Analyses  of  a  Monolith 
Using  SAP2000  Linear  Elastic 
Finite  Elements 


A  dynamic  time-history  structural  analysis  was  performed  of  the  generator 
bay  and  erection/service  bay  substructure  finite  element  models  in  order  to 
calculate  the  top-of-substructure  acceleration  time-history  and  then  the  pseudo- 
spectral  acceleration  curve.  The  analytical  model  and  the  dynamic  earthquake 
analyses  were  performed  using  the  structural  analysis  software  SAP2000 
(Computers  and  Structures,  Inc.,  19971).  The  analytical  model  represents  the 
geometry  and  material  properties  of  the  idealized  powerhouse  substructures. 
Parametric  analyses  were  carried  out  to  determine  the  type  of  finite  elements  to 
be  used  in  the  study. 

The  response  of  elements  to  applied  loads  or  displacements  will  depend  on 
the  material  constitutive  properties.  The  constitutive  properties  relate  the  stresses 
to  strains  as  shown  in  Equation  C-l : 

M=[c]{s)  (C-l) 


where  {cr}  is  the  stresses  vector,  {f}  is  the  strains  vector,  and  [c]  is  the  con¬ 
stitutive  matrix.  Equations  C-2  and  C-3  list  the  constitutive  matrix  for  two- 
dimensional  (2-D)  plane  stress  elements  and  for  2-D  plane  strain  elements, 
respectively.  For  linear  elastic  structural  analysis,  it  is  typically  assumed  that  the 
material  constitutive  matrix  [c]  doesn’t  change  with  load. 


1  References  cited  in  this  appendix  are  included  in  the  References  section  at  the  end  of  the  main 
text. 
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Cl 


1 


V 


0 

0 


(C-2) 


M-l 


-  ir 


u  1 


0  0 


(l-o) 


[C]  = 


(l  +  L>Xl  -  2 u) 


l-v  v  0 

o  l-o  0 

0  0  UN 


where  E  is  the  elastic  modulus  and  o  is  Poisson’s  ratio. 


(C-3) 


Two  types  of  plane  elements  are  available  in  SAP2000  to  conduct  a  planar 
analysis:  plane  stress  elements  or  plane  strain  elements.  When  plane  stress 
elements  are  used,  it  is  assumed  that  the  stresses  out  of  plane  to  the  element  are 
zero  (gz  =  0);  but  out-of-plane  deformations  will  occur.  With  the  plane  strain 
elements,  the  out-of-plane  deformations  are  zero  (sz  =  0),  and  out-of-plane 
stresses  are  nonzero. 


Monoliths  of  concrete  gravity  dams  are  constructed  with  construction  joints 
between  them,  allowing  out-of-plane  deformations.  Plane  stress  elements  are 
recommended  in  this  situation  (HQUSACE  2003).  A  similar  construction  process 
is  applied  to  powerhouse  monoliths,  so  it  is  logical  to  make  a  similar  modeling 
assumption. 

In  SAP2000  structural  analysis  software,  different  numbers  of  nodes  can  be 
assigned  to  planar  finite  elements.  The  number  of  nodes  will  define  the  total 
degrees  of  freedom  for  an  element.  In  this  parametric  study,  plane  stress  and 
plane  strain  elements  with  four  and  nine  nodes  and  shell  elements  were  com¬ 
pared.  The  1 1 0-ft-height  erection/service  bay  block  model  was  used.  The  first 
mode  period,  the  number  of  modes  that  must  be  included  based  on  participating 
mass  ratio  (PMR)  considerations,  and  the  peak  acceleration  at  top  of  the  sub¬ 
structure  were  all  compared.  The  investigation  was  performed  using  the  central- 
eastern  United  States  (CEUS)  and  the  Pacific  Northwest  of  the  United  States 
(PNW)  acceleration  time-histories. 

Table  C-l  lists  the  first  mode  periods  and  frequencies,  highest  periods  and 
frequencies,  the  number  of  modes  included,  and  the  resulting  PMR  for  the 
different  elements  considered  in  this  study  for  the  1 1 0-ft-height  erection/service 
bay  block  model.  Significant  differences  in  the  model  dynamic  properties  were 
not  found  between  the  types  of  elements.  The  first  mode  period  for  the  stress  and 
shell  elements  were  0.091  sec  and  for  the  strain  elements  0.090  sec  for  a  percent 
difference  of  1.09.  There  was  no  significant  difference  in  the  first  mode  periods 
for  the  four-node  and  nine-node  elements  for  the  same  type  of  planar  element 
(i.e.,  stress  or  strain).  Tables  C-2  and  C-3  list  the  top-of-substructure  total  and 
relative  peak  acceleration  for  the  CEUS  and  PNW  ground  motions,  respectively. 
No  significant  differences  were  found  between  the  plane  strain  and  plane  stress 
elements  given  the  same  number  of  nodes.  In  addition,  no  significant  difference 
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was  observed  between  plane  stress  (or  plane  strain)  when  the  number  of  nodes 
was  increased  for  each  element.  Note  that  these  observations  will  be  dependent 
on  the  discretization  of  the  fine  mesh,  as  reflected  by  a  large  number  of  vertical 
layers  of  elements  in  the  mesh. 


Table  C-1 

Periods  and  Frequencies  Found  for  Different  Types  of  Elements  for  110-ft-Height 
Erection/Service  Substructure  Model 

Model  ID 

Vertical  No. 
of  Layers 
of 

Elements  in 
Mesh 

Element 

Type 

Modes  / 

Period,  sec 

Frequency,  Hz  I 

T, 

Ti 

fi 

f. 

4ShErD110B5 

4 

Shell 

4 

93.2 

0.0917 

0.0298 

10.91 

33.47 

8ShErD110B5 

8 

Shell 

6 

92.2 

0.0914 

0.0216 

10.94 

46.22 

4S4ErD110B5 

4 

Q4  Stress 

4 

93.3 

0.0910 

0.0289 

10.98 

34.6 

8S4ErD110B5 

8 

Q4  Stress 

6 

92.4 

0.0912 

0.0216 

10.96 

46.30 

16S4ErD110B5 

16 

Q4  Stress 

6 

90.3 

0.0914 

0.0215 

10.94 

46.59 

4S9ErD110B5 

4 

Q9  Stress 

6 

93.2 

0.0914 

0.0214 

10.94 

46.75 

8S9ErD110B5 

8 

Q9  Stress 

8 

93.3 

0.0915 

0.0155 

10.93 

64.43 

4N4ErD110B5 

4 

Q4  Strain 

4 

93.4 

0.0898 

0.0284 

11.13 

35.22 

8N4ErD110B5 

8 

Q4  Strain 

6 

92.5 

0.0901 

0.0215 

11.10 

46.49 

16N4ErD110B5 

16 

Q4  Strain 

6 

90.5 

0.0903 

0.0214 

11.08 

46.77 

4N9ErD110B5 

4 

Q9  Strain 

6 

91.8 

0.0903 

0.0213 

11.08 

46.77 

8N9ErD110B5 

8 

Q9  Strain 

8 

93.4 

0.0904 

0.0153 

11.07 

65.12 
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Table  C-2 

Absolute  and  Relative  Peak  Acceleration  at  Top-of-substructure  for  110-ft-Height 
Erection/Service  Bay  Block  Model  Using  Different  Types  of  Finite  Elements  for  CEUS 
Ground  Motion 

Model  ID 

Vertical  No.  of 
Layers  of 
Elements  in 

Mesh 

Element  Type 

Total  Peak  Acceleration 

9 

Relative  Peak  Acceleration 

g 

(+) 

(■) 

(+) 

(-) 

4ShErD110B5 

4 

Shell 

0.44 

-0.54 

0.44 

-0.52 

8ShErD110B5 

8 

Shell 

0.46 

-0.56 

0.46 

-0.53 

4S4ErD110B5 

4 

Q4  Stress 

0.44 

-0.54 

0.44 

-0.51 

8S4ErD110B5 

8 

Q4  Stress 

0.46 

-0.55 

0.45 

-0.53 

16S4ErD110B5 

16 

Q4  Stress 

0.46 

-0.56 

0.46 

-0.54 

4S9ErD110B5 

4 

Q9  Stress 

0.46 

-0.56 

0.46 

-0.54 

8S9ErD110B5 

8 

Q9  Stress 

0.46 

-0.56 

0.46 

-0.54 

4N4ErD110B5 

4 

Q4  Strain 

0.44 

-0.52 

0.45 

-0.50 

8N4ErD110B5 

8 

Q4  Strain 

0.46 

-0.54 

0.46 

-0.52 

16N4ErD110B5 

16 

Q4  Strain 

0.46 

-0.55 

0.47 

-0.53 

4N9ErD110B5 

4 

Q9  Strain 

0.46 

-0.55 

0.47 

-0.53 

8N9ErD110B5 

8 

Q9  Strain 

0.46 

-0.55 

0.46 

-0.54 

Table  C-3 

Absolute  and  Relative  Peak  Acceleration  at  Top-of-substructure  for  110-ft-Height 
Erection/Service  Bay  Block  Model  Using  Different  Types  of  Finite  Elements  for  PNW 
Ground  Motion 

Model  ID 

Vertical  No.  of 
Layers  of 
Elements  in 

Mesh 

Element  Type 

Total  Peak  Acceleration 

g 

Relative  Peak  Acceleration 

g 

(+) 

(■) 

(+) 

(-) 

4ShErD110B5 

4 

Shell 

0.40 

-0.35 

0.34 

-0.33 

8ShErD110B5 

8 

Shell 

0.41 

-0.36 

0.36 

-0.35 

4S4ErD110B5 

4 

Q4  Stress 

0.40 

-0.35 

0.35 

-0.34 

8S4ErD110B5 

8 

Q4  Stress 

0.41 

-0.36 

0.36 

-0.35 

16S4ErD110B5 

16 

Q4  Stress 

0.41 

-0.37 

0.36 

-0.35 

4S9ErD110B5 

4 

Q9  Stress 

0.41 

-0.37 

0.36 

0.35 

8S9ErD110B5 

8 

Q9  Stress 

0.41 

-0.37 

0.36 

0.35 

4N4ErD110B5 

4 

Q4  Strain 

0.41 

-0.37 

0.35 

-0.35 

8N4ErD110B5 

8 

Q4  Strain 

0.41 

-0.37 

0.36 

-0.36 

16N4ErD110B5 

16 

Q4  Strain 

0.42 

-0.37 

0.36 

0.36 

4N9ErD110B5 

4 

Q9  Strain 

0.41 

-0.37 

0.36 

-0.36 

8N9ErD110B5 

8 

Q9  Strain 

0.42 

0.37 

0.36 

-0.35 
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Appendix  D 

Hydrodynamic  Added  Mass 


D.l.  Introduction 

Hydrodynamic  added  mass  was  added  to  the  upstream  and  downstream  faces 
of  the  finite  element  models  to  account  for  the  pool  and  tailrace  water  elevations, 
respectively,  of  the  powerhouse  monoliths.  Specifically,  the  mass  was  added  to 
the  upstream  and  downstream  nodes.  For  the  computations  of  the  water-added 
mass,  a  simplified  procedure  of  the  Westergaard’s  Solution  (1993) 1  found  in 
HQUSACE  (1995  a)  was  adapted  for  this  study. 

The  simplified  procedure  assumes  a  parabolic  mass  distribution  along  the 
upstream  and  downstream  faces  of  the  powerhouse  substructure  model  as  shown 
in  Figure  D-l.  In  order  to  calculate  the  added  mass,  the  effective  width  of  water 
parallel  to  the  monolith  at  the  point  of  interest  at  level  (y)  is  calculated.  The 
effective  width  can  be  approximated  by  Equation  D-l: 

=  (D-l) 

O 


where  Hw  is  the  total  depth  of  water  and  y  the  distance  below  the  free  surface  of 
the  water  to  the  point  of  interest.  This  “width”  times  the  surface  area  associated 
with  the  node  estimates  the  volume  of  water  and  consequently  the  mass  to  be 
added  to  the  node.  This  calculation  is  repeated  for  all  nodes  along  the  upstream 
and  downstream  faces  of  the  substructure.  The  total  added  weight  is  calculated 
by  Equation  D-2: 

Wtotal-^wwHw2  (D-2) 

per  Ebeling  and  Morrison  (1992)  where  yw  is  the  unit  weight  of  water  and  w  is 
the  plan  width  (i.e.,  along  the  axis  of  the  powerhouse  substructure)  of  the 
idealized  substructure  monolith.  The  weight  is  then  divided  by  the  gravitational 
constant  of  acceleration  g  to  get  the  total  added  mass  for  the  model. 


1  References  cited  in  this  appendix  are  included  in  the  References  section  at  the  end  of  the  main 
text. 
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Figure  D-1.  Hydrodynamic  added  mass  distribution  on  the  idealized  powerhouse 
substructure  model 

D.2.  Example  Computation 

The  following  is  an  example  of  the  calculations  for  the  1 1 0-ft-height  block 
model  of  the  erection/service  bay  idealized  substructure  (Figure  D-2).  The  finite 
element  model  has  a  height  of  1 10  ft  and  includes  33  nodes  along  the  upstream 
face  and  23  nodes  along  the  downstream  face.  The  nodes  are  spaced  at  intervals 
of  3.4375  ft.  The  model  has  a  plan  width  of  75  ft  (along  the  axis  of  the  power¬ 
house  substructure),  the  upstream  pool  water  elevation  is  assigned  20  ft  above  the 
top  of  the  idealized  substructure,  and  the  tailrace  water  elevation  is  34.4  ft  below 
the  top  of  the  idealized  substructure  model.  The  unit  weight  of  water  was  taken 
as  62.4  lb/ft3. 

Table  D-1  summarizes  the  calculations  for  the  pool  added  mass  in  spread¬ 
sheet  form.  In  the  first  column  the  node  numbers  are  identified  and  in  Column  2 
the  distance  below  the  water  table  of  each  node  is  calculated.  The  effective  width 
b(y)  is  calculated  in  Column  3  using  Equation  D- 1 .  The  added  weight  at  each 
node  is  calculated  by  multiplying  the  effective  width  b(y)  by  the  tributary  area  of 
each  node  times  the  water  unit  weight  as  presented  in  Column  4.  For  the  node  at 
top  of  the  model,  the  20  ft  of  water  weight  above  the  top-of-substructure  was 
added.  The  weight  is  divided  by  the  gravitational  acceleration  constant 
32.2  fit/sec2  in  Column  5  to  obtain  the  added  hydrodynamic  mass  at  each  node. 
Table  D-2  lists  the  results  for  the  tailrace  water  added  mass  for  each  downstream 
face  node. 

The  total  weight  is  calculated  by  summing  the  weight  for  each  of  the  nodes. 
This  value  is  compared  to  the  exact  value  calculated  by  Equation  D-2.  Accept¬ 
able  errors  of  1.52  percent  were  found  for  the  pool  side  and  1.02  percent  for  the 
tailrace  side  (this  error  results  from  the  linearization  of  the  Figure  D-1  curve 
between  nodes). 
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Figure  D-2.  Finite  element  block  model  for  1 1 0-ft-height  erection/service  bay 
idealization 

Table  D-3  lists  the  percent  contribution  of  the  added  mass  to  the  model  mass. 
For  this  example  28.2  percent  of  the  total  mass  in  the  block  finite  element  model 
is  attributed  to  the  hydrodynamic  added  mass. 
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Table  D-1 

Calculated  Hydrodynamic  Mass  Added  to  Each  Node  along  the 
Upstream  Pool  Face 

i 

2 

3 

4 

5 

y 

b 

Weight 

Mass 

Node 

ft 

ft 

lb 

(lb-sec2)/ft 

33 

20 

44.61642074 

2446932 

75991.67159 

32 

23.4375 

48.29869224 

777005.2 

24130.59663 

31 

26.875 

51.71945808 

832036.8 

25839.65161 

30 

30.3125 

54.92759794 

883647.7 

27442.47615 

29 

33.75 

57.95843069 

932406.3 

28956.71595 

28 

37.1875 

60.83846119 

978738.7 

30395.61318 

27 

40.625 

63.58818311 

1022975 

31769.4067 

26 

44.0625 

66.22383033 

1065376 

33086.20715 

25 

47.5 

68.7585222 

1106153 

34352.56913 

24 

50.9375 

71.20304099 

1145479 

35573.87957 

23 

54.375 

73.56637628 

1183499 

36754.62976 

22 

57.8125 

75.85611637 

1220335 

37898.61093 

21 

61.25 

78.0787363 

1256092 

39009.05808 

20 

64.6875 

80.2398139 

1290858 

40088.75795 

19 

68.125 

82.34419473 

1324712 

41140.13145 

18 

71.5625 

84.39612004 

1357723 

42165.29755 

17 

75 

86.39932725 

1389949 

43166.12352 

16 

78.4375 

88.35713 

1421445 

44144.26487 

15 

81.875 

90.27248248 

1452259 

45101.19758 

14 

85.3125 

92.14803181 

1482431 

46038.24415 

13 

88.75 

93.98616088 

1512002 

46956.59513 

12 

92.1875 

95.78902395 

1541006 

47857.3268 

11 

95.625 

97.55857615 

1569474 

48741.41596 

10 

99.0625 

99.2965984 

1597434 

49609.75238 

9 

102.5 

101.0047183 

1624913 

50463.14926 

8 

105.9375 

102.6844282 

1651936 

51302.35214 

7 

109.375 

104.3371002 

1678523 

52128.04654 

6 

112.8125 

105.9639993 

1704696 

52940.86454 

5 

116.25 

107.566295 

1730473 

53741.39038 

4 

119.6875 

109.1450708 

1755871 

54530.16543 

3 

123.125 

110.7013331 

1780908 

55307.69244 

2 

126.5625 

112.2360184 

1805597 

56074.43932 

1 

130 

113.75 

914976.6 

28415.4212 

Total 

45437862 

1411113.715 

Exact  Value 

46137000 

1432826.087 

%  Error  = 

1.52 

1.52 
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Table  D-2 

Calculated  Hydrodynamic  Mass  Added  to  Each  Node  along  the 
Downstream  Tailrace  Face 

i 

2 

3 

4 

5 

y 

b 

Weight 

Mass 

Node 

ft 

ft 

lb 

(lb-sec2)/ft 

1079 

0 

0 

0 

0 

1078 

3.436364 

14.10322739 

113405.3 

3521.904768 

1077 

6.872727 

19.94497544 

320758.7 

9961.450977 

1076 

10.30909 

24.42750638 

392847.6 

12200.236 

1075 

13.74545 

28.20645477 

453621.3 

14087.61907 

1074 

17.18182 

31.53577514 

507164.1 

15750.43694 

1073 

20.61818 

34.54571082 

555570.4 

17253.73921 

1072 

24.05455 

37.31363235 

600084.6 

18636.16831 

1071 

27.49091 

39.88995089 

641517.4 

19922.90195 

1070 

30.92727 

42.30968216 

680432 

21131.42861 

1069 

34.36364 

44.5983209 

717238.3 

22274.48154 

1068 

37.8 

46.77511358 

752245.9 

23361.67333 

1067 

41.23636 

48.85501277 

785695.2 

24400.47199 

1066 

44.67273 

50.84990949 

817777.5 

25396.81646 

1065 

48.10909 

52.76944493 

848647.8 

26355.52198 

1064 

51.54545 

54.6215648 

878433.9 

27280.55703 

1063 

54.98182 

56.41290955 

907242.7 

28175.23814 

1062 

58.41818 

58.14909618 

935164.3 

29042.37072 

1061 

61.85455 

59.83492633 

962276.2 

29884.35293 

1060 

65.29091 

61.47454296 

988644.8 

30703.25395 

1059 

68.72727 

63.07155028 

1014328 

31500.87389 

1058 

72.16364 

64.62910704 

1039377 

32278.79038 

1057 

75.6 

66.15 

531918.2 

16519.19763 

Total 

15444391 

479639.4858 

Exact  Value 

15602933 

484563.1304 

%  Error  = 

1.02 

1.02 

Table  D-3 

Summary  of  Hydrodynamic  Mass  Added  to  the  Idealized  Model 

Hydrodynamic 
Added  Mass 

Block 

(k-sec2)/ft 

Pool 

(k-sec2)/ft 

Tail 

(k-sec2)/ft 

Total  Added 

Mass 

(k-sec2)/ft 

Total 

(k-sec2)/ft 

Mass 

4803.97 

1411.1 

479.6 

1890.7 

6694.67 

%  of  Added  Mass 

74.6% 

25.4% 

100.0% 

%  of  Block  Mass 

29.4% 

10.0% 

39.4% 

%  of  Total  Mass 

21.1% 

7.2% 

28.2% 
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Appendix  E 

Absolute  Acceleration  versus 
Pseudo-Spectral  Acceleration 
(PSA) 


This  appendix  compares  amplification  factor  curves  developed  using  abso¬ 
lute  acceleration  (SA)  curves  and  pseudo-spectral  acceleration  (PSA)  curves.  The 
1 1 0-ft-height  block  model  of  the  erection/service  bay  substructure  is  used  in  this 
example. 

A  response  spectrum  is  a  graph  of  maximum  values  of  acceleration,  velocity, 
or  displacement  response  of  an  infinite  series  of  damped  elastic  single-degree-of- 
freedom  (SDOF)  systems  subjected  to  an  acceleration  time-history  (Ebeling 
1992). 1  When  referring  to  absolute  acceleration,  SA  (also  known  as  spectral 
acceleration  response  spectrum),  the  ground  acceleration  time-history  is  added  to 
the  relative  response  values  of  the  computed  time-history  (refer  to  SA  equation 
given  in  box  in  Figure  1-3).  For  the  PSA,  the  computed  relative  time -history 
acceleration  is  used  to  construct  the  response  spectrum  curve. 

Figure  E-l  shows  the  calculated  top-of-substructure  response  spectra  with 
5  percent  damping  using  SA  and  PSA  for  the  1 1 0-ft-height  block  model.  The 
structural  analysis  software  SAP2000  (Computers  and  Structures,  Inc.,  1997)  was 
used  for  the  calculation  of  the  response  spectra,  and  the  central-eastern  United 
States  ground  motion  (CEUS)  as  the  earthquake.  Figure  E-2  shows  the  resulting 
amplification  factor  curves  constructed  using  the  SA  and  the  PSA  curves  for 
5  percent  damping. 

The  difference  between  the  two  spectra  and  amplification  factor  curves  is 
small  because  a  short-period  structure  with  low  percent  damping  was  being 
evaluated.  The  PSA  considers  only  elastic  forces  and  SA  considers  elastic  and 
damping  forces  (Appendix  A,  Ebeling  1992).  With  low  percent  damping,  the 
contribution  of  damping  forces  considered  in  SA  is  not  significant  to  the  analysis. 
For  long-period  structures  and  a  high  percent  damping,  the  difference  between 
the  two  spectra  will  be  increased. 


1  References  cited  in  this  appendix  are  included  in  the  References  section  at  the  end  of  the  main 
text. 
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Figure  E-1 .  SA  and  PSA  response  spectrum  curves  at  top-of-substructure  for  the 
1 1 0-ft-height  block  model  using  the  CEUS  ground  motion  for 
5  percent  damping 


Figure  E-2.  AF  using  SA  and  PSA  at  top-of-substructure  for  the  1 1 0-ft-height 
block  model  using  the  CEUS  ground  motion 
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Appendix  F 

Amplification  Factor  Curves 
for  2  Percent  Damping 


F.l.  Introduction 

In  this  appendix,  amplification  factor  (AF)  curves  for  2  percent  damping  are 
shown  for  the  composite  model  of  the  idealized  generator  bay  substructure.  The 
results  for  the  block  model  and  block- frame-shear  wall  model  of  the  erection/ 
service  bay  substructure  are  also  presented.  Lower  values  of  damping  are  used  to 
compute  response  spectra  for  analyzing  the  equipment  on  top  of  the  substructure 
than  are  used  for  analyzing  the  superstructure. 

The  AF  curves  were  calculated  using  the  same  procedure  as  explained  in 
Chapter  1 .  The  superstructure  was  decoupled  from  the  substructure,  and  a  time- 
history  analysis  was  performed  of  the  four  idealized  substructure  finite  element 
models.  The  top-of-substructure  2  percent  damping  response  spectrum  was  then 
calculated.  The  percent  damping  for  the  substructure  finite  element  model  was 
kept  at  5  percent  for  all  analyses.  To  calculate  the  AF  curve,  the  top-of- 
substructure  spectral  acceleration  for  2  percent  damping  was  divided  by  the 
pseudo-spectral  acceleration  in  rock  for  2  percent  damping  for  each  frequency. 
The  time-histories  of  the  ground  motions  for  the  central-eastern  United  States 
(CEUS)  and  the  Pacific  Northwest  of  the  United  States  (PNW)  were  used  in  the 
analyses.  The  AF  curves  presented  in  this  appendix  are  plotted  up  to  50  Hz. 
However,  the  ground  motions  used  in  this  study  have  frequency  content  up  to 
33  Hz.  The  use  of  AF  curves  above  33  Hz  is  not  recommended. 

Figure  F-l  shows  the  process  used  to  develop  the  AF  curve  and  its  applica¬ 
tion  when  doing  response  spectrum  analysis  for  equipment  resting  on  the  power¬ 
house  substructure. 
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a.  Top-of-substructure  2  percent  damping  response  spectrum  calculation 
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b.  Construction  of  the  AF  curve  for  2  percent  damping 

Figure  F-1 .  AF  curve  development  for  2  percent  damping  and  its  application  to  equipment  analysis 
(Continued) 
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F.2.  Generator  Bay  Substructure  Amplification 
Factor  for  2  Percent  Damping 

Table  F-l  lists  the  AF  curves  and  response  spectrum  figures  obtained  for  the 
composite  structural  model  of  the  idealized  generator  bay  substructure.  As 
expected,  when  the  percent  damping  value  is  decreased,  the  response  increases 
for  each  of  the  four  group  models.  Tables  F-2  and  F-3  list  the  SA  and  PGA  ratio 
at  50  Hz  for  dry  and  wet  conditions,  respectively. 


Table  F-1 

List  of  Figures  for  Generator  Bay  Substructure  AF  Curves 

Group 

Height,  ft 

CEUS 

PNW 

1 

125 

Figure  F-2 

Figure  F-6 

2 

100 

Figure  F-3 

Figure  F-7 

3 

70 

Figure  F-4 

Figure  F-8 

4 

40 

Figure  F-5 

Figure  F-9 
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Figure  F-2.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
125-ft  height  (Dry  Model  ID:  25S9GeD125C5,  Wet  Model  ID: 
25S9GeW125C5),  for  CEUS  ground  motion  with  2  percent  damping 
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a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


b.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 
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c.  CEUS  ground  motion  response  spectrum 


Figure  F-3.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
100-ft  height  (Dry  Model  ID:  50S9GeD100C5,  Wet  Model  ID: 
50S9GeW100C5),  for  CEUS  ground  motion  with  2  percent  damping 
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Figure  F-4.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
70-ft  height  (Dry  Model  ID:  35S9GeD70C5,  Wet  Model  ID: 
35S9GeW70C5),  for  CEUS  ground  motion  with  2  percent  damping 
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b.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 
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Figure  F-5.  Response  spectra  and  AF  curves,  generator  bay,  composite  model 
40-ft  height  (Dry  Model  ID:  8S9GeD40C5,  Wet  Model  ID: 
8S9GeW40C5),  for  CEUS  ground  motion  with  2  percent  damping 
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a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


Figure  F-6.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
125-ft  height  (Dry  Model  ID:  25S9GeD125C5,  Wet  Model  ID: 
25S9GeW125C5),  for  PNW  ground  motion  with  2  percent  damping 
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Figure  F-7.  Response  spectra  and  AF  curves,  generator  bay,  composite  model 
100-ft  height  (Dry  Model  ID:  50S9GeD100C5,  Wet  Model  ID: 
50S9GeW100C5),  for  PNW  ground  motion  with  2  percent  damping 
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c.  PNW  ground  motion  response  spectrum 

Figure  F-8.  Response  spectra  and  AF  curves,  generator  bay,  composite  model, 
70-ft  height  (Dry  Model  ID:  35S9GeD70C5,  Wet  Model  ID: 
35S9GeW70C5),  for  PNW  ground  motion  with  2  percent  damping 
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Figure  F-9.  Response  spectra  and  AF  curves,  generator  bay,  composite  model 
40-ft  height  (Dry  Model  ID:  8S9GeD40C5,  Wet  Model  ID: 
8S9GeW40C5),  for  PNW  ground  motion  with  2  percent  damping 
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Table  F-2 

SA  at  Frequency  of  50  Hz  for  Generator  Bay  Composite  Models  for 

2  Percent  Damping,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g 

sa5„Hz/pga 

g 

SA50HZ 

g 

SA5ohz/PGA 

g 

25S9GeD125C5 

125 

0.62 

3.26 

0.44 

2.75 

50S9GeD100C5 

100 

0.51 

2.68 

0.30 

1.88 

35S9GeD70C5 

70 

0.4 

2.11 

0.23 

1.44 

8S9GeD40C5 

40 

0.3 

1.58 

0.18 

1.13 

Table  F-3 

SA  at  Frequency  of  50  Hz  for  Generator  Bay  Composite  Models  for 

2  Percent  Damping,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g 

SA50Hz/PGA 

g 

SA50HZ 

g 

SAsohz/PGA 

g 

25S9GeW125C5 

125 

0.54 

2.84 

0.42 

2.63 

50S9GeW100C5 

100 

0.55 

2.89 

0.38 

2.38 

35S9GeW70C5 

70 

0.37 

1.95 

0.29 

1.81 

8S9GeW40C5 

40 

0.26 

1.37 

0.19 

1.1875 

F.3.  Block  Model  of  Idealized  Erection/Service 
Bay  Substructure  AF  Curves  for  2  Percent 
Damping 

The  AF  curves  for  the  block  model  with  5  percent  damping  were  shown  in 
Chapter  3.  In  this  section,  the  AF  and  response  spectrum  curves  are  shown  for 
2  percent  damping.  As  for  5  percent  damping,  wet  and  dry  conditions  were 
considered  and  two  earthquakes  were  used  to  analyze  the  idealized  models: 
CEUS  ground  motion  and  PNW  ground  motion.  Table  F-4  lists  the  figures  for 
each  group  and  ground  motion  used  in  the  analysis.  The  AFs  obtained  with 
2  percent  damping  were  higher  than  the  results  with  5  percent  damping,  as 
expected.  Tables  F-5  and  F-6  list  that  SA  and  PGA  ratio  for  dry  and  wet 
conditions,  respectively. 


Table  F-4 

List  of  Figures  for  Substructure  Block  Model  AF  Curves 

Group 

Height 

ft 

CEUS 

PNW 

1 

110 

Figure  F-10 

Figure  F-14 

2 

75 

Figure  F-11 

Figure  F-15 

3 

45 

Figure  F-12 

Figure  F-16 

4 

20 

Figure  F-13 

Figure  F-17 
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Figure  F-10.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  110-ft  height  (Dry  Model  ID:  16S9ErD110B5,  Wet  Model  ID: 
16S9ErW110B5),  for  CEUS  ground  motion  with  2  percent  damping 
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.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 
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Figure  F-11.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  75-ft  height  (Dry  Model  ID:  16S9ErD75B5,  Wet  Model  ID: 
16S9ErW75B5),  for  CEUS  ground  motion  with  2  percent  damping 
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Figure  F-12.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  45-ft  height  (Dry  Model  ID:  16S9ErD45B5,  Wet  Model  ID: 
16S9ErW45B5),  for  CEUS  ground  motion  with  2  percent  damping 
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Figure  F-13.  Response  spectra  and  AF  curves,  erection/service  bay,  block 
model,  20-ft  height  (Dry  Model  ID:  8S9ErD20B5,  Wet  Model  ID: 
8S9ErW20B5),  for  CEUS  ground  motion  with  2  percent  damping 
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c.  PNW  ground  motion  response  spectrum 


Figure  F-14.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  110-ft  height  (Dry  Model  ID:  16S9ErD110B5,  Wet  Model  ID: 
16S9ErW110B5),  for  PNW  ground  motion  with  2  percent  damping 
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Figure  F-15.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  75-ft  height  (Dry  Model  ID:  16S9ErD75B5,  Wet  Model  ID: 
16S9ErW75B5),  for  PNW  ground  motion  with  2  percent  damping 
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Figure  F-16.  Response  spectra  and  AF  curves,  erection/service  bay,  block 

model,  45-ft  height  (Dry  Model  ID:  16S9ErD45B5,  Wet  Model  ID: 
16S9ErW45B5),  for  PNW  ground  motion  with  2  percent  damping 
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Figure  F-17.  Response  spectra  and  AF  curves,  erection/service  bay,  block 
model,  20-ft  height  (Dry  Model  ID:  8S9ErD20B5,  Wet  Model  ID: 
8S9ErW20B5),  for  PNW  ground  motion  with  2  percent  damping 
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Table  F-5 

SA  at  Frequency  of  50  Hz  for  Generator  Bay  Composite  Models  for 

2  Percent  Damping,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g 

sa50Hz/pga 

g 

SA50HZ 

g 

SA50Hz/PGA 

g 

16S9ErD110B5 

110 

0.58 

3.05 

0.41 

2.56 

16S9ErD75B5 

75 

0.4 

2.11 

0.23 

1.44 

16S9ErD45B5 

45 

0.28 

1.47 

0.17 

1.06 

8S9ErD20B5 

20 

0.24 

1.26 

0.16 

1.00 

Table  F-6 

SA  at  Frequency  of  50  Hz  for  Generator  Bay  Composite  Models  for 

2  Percent  Damping,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g 

sa50Hz/pga 

g 

SA50HZ 

g 

SAsohz/PGA 

g 

16S9ErW110B5 

110 

0.68 

3.58 

0.44 

2.75 

16S9ErW75B5 

75 

0.4 

2.11 

0.31 

1.94 

16S9ErW45B5 

45 

0.26 

1.37 

0.17 

1.06 

8S9ErW20B5 

20 

0.2 

1.05 

0.16 

1.00 

F.4.  Block-Frame-Shear  Wall  Model  of  Idealized 
Erection/Service  Bay  Substructure  Amplification 
Factor  Curves  for  2  Percent  Damping 

In  Section  F.2,  the  AF  curves  were  shown  for  the  composite  structural  model 
of  the  generator  bay  substructure.  For  the  erection/service  bay  substructure  two 
models  were  analyzed:  block  model  and  block-frame-shear  wall  model.  The  AF 
curves  for  the  block  model  were  shown  in  the  previous  section.  In  this  section  the 
AF  and  response  spectrum  curves  for  2  percent  damping  are  presented  for  the 
block- frame-shear  wall  model.  Table  F-7  lists  the  figures  shown  in  this  section 
for  the  two  ground  motions  considered  in  this  section.  The  models  were  analyzed 
assuming  dry  and  wet  conditions.  Tables  F-8  and  F-9  list  the  SA  and  PGA  ratios 
at  50  Hz  for  dry  and  wet  conditions,  respectively. 


Table  F-7 

List  of  Figures  for  Substructure  Block  Model  AF  Curves 

Group 

Height 

ft 

CEUS 

PNW 

1 

110 

Figure  F-18 

Figure  F-22 

2 

75 

Figure  F-19 

Figure  F-23 

3 

45 

Figure  F-20 

Figure  F-24 

4 

20 

Figure  F-21 

Figure  F-25 
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Figure  F-18.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  110-ft  height  (Dry  Model  ID:  14S9ErD110FW5, 
Wet  Model  ID:  14S9ErW1 10FW5),  for  CEUS  ground  motion  with 
2  percent  damping 
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Figure  F-19.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame 
shear  wall  model,  75-ft  height  (Dry  Model  ID:  20S9ErD75FW5,  Wet 
Model  ID:  20S9ErW75FW5),  for  CEUS  ground  motion  with 
2  percent  damping 


Appendix  F  Amplification  Factor  Curves  for  2  Percent  Damping 


SA(g)  q-  SA(g)  (Amplification  Factor)SA 


Frequency  (Hz) 

a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 


0  5  10  15  20  25  30  35  40  45  50 

Frequency  (Hz) 

.  Response  spectrum  at  top-of-substructure  for  wet  and  dry  conditions 


0  5  10  15  20  25  30  35  40  45  50 

Frequency  (Hz) 
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Figure  F-20.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  45-ft  height  (Dry  Model  ID:  12S9ErD45FW5,  Wet 
Model  ID:  12S9ErW45FW5),  for  CEUS  ground  motion  with 
2  percent  damping 
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Figure  F-21.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  20-ft  height  (Dry  Model  ID:  8S9ErD20FW5,  Wet 
Model  ID:  8S9ErW20FW5),  for  CEUS  ground  motion  with  2  percent 
damping 
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Figure  F-22.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  110-ft  height  (Dry  Model  ID:  14S9ErD110FW5, 
Wet  Model  ID:  14S9ErW1 10FW5),  for  PNW  ground  motion  with 
2  percent  damping 
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Figure  F-23.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  75-ft  height  (Dry  Model  ID:  20S9ErD75FW5,  Wet 
Model  ID:  20S9ErW75FW5),  for  PNW  ground  motion  with  2  percent 
damping 
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a.  AF  at  top-of-substructure  for  wet  and  dry  conditions 
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Figure  F-24.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  45-ft  height  (Dry  Model  ID:  12S9ErD45FW5,  Wet 
Model  ID:  12S9ErW45FW5),  for  PNW  ground  motion  with  2  percent 
damping 
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Figure  F-25.  Response  spectra  and  AF  curves,  erection/service  bay,  block-frame- 
shear  wall  model,  20-ft  height  (Dry  Model  ID:  8S9ErD20FW5,  Wet 
Model  ID:  8S9ErW20FW5),  for  PNW  ground  motion  with  2  percent 
damping 


Appendix  F  Amplification  Factor  Curves  for  2  Percent  Damping 


Table  F-8 

SA  at  Frequency  of  50  Hz  for  Generator  Bay  Composite  Models  for 

2  Percent  Damping,  Dry  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50HZ 

g 

SAsohz/PGA 

g 

SA50HZ 

g 

SA5ohz/PGA 

g 

14S9ErD110FW5 

110 

0.63 

3.32 

0.41 

2.56 

20S9ErD75FW5 

75 

0.41 

2.16 

0.31 

1.94 

12S9ErD45FW5 

45 

0.36 

1.89 

0.18 

1.13 

8S9ErD20FW5 

20 

0.22 

1.16 

0.17 

1.0625 

Table  F-9 

SA  at  Frequency  of  50  Hz  for  Generator  Bay  Composite  Models  for 

2  Percent  Damping,  Wet  Condition 

Model  ID 

Height 

ft 

CEUS 

PNW 

SA50Hz 

g 

sa50Hz/pga 

g 

SA50Hz 

g 

SA50Hz/PGA 

g 

14S9ErW110FW5 

110 

0.6 

3.16 

0.44 

2.75 

20S9ErW75FW5 

75 

0.47 

2.47 

0.32 

2.00 

12S9ErW45FW5 

45 

0.3 

1.58 

0.21 

1.31 

8S9ErD20FW5 

20 

0.33 

1.74 

0.17 

1.0625 

F30 


Appendix  F 


Amplification  Factor  Curves  for  2  Percent  Damping 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington, 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not 
display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. _ 


1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

June  2006  Final  report 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Response  Amplification  of  Idealized  Powerhouse  Substructures  to  Earthquake  Ground 
Motions 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Robert  M.  Ebeling,  Elvin  J.  Perez-Marcial,  Donald  E.  Yule 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center 

Information  Technology  Laboratory  and  Geotechnical  and  Structures  Laboratory 
3909  Halls  Ferry  Road,  Vicksburg,  MS  39180-6199; 

CMA  Architects  and  Engineers,  LLP 

1509  F.D.  Roosevelt  Avenue,  Guaynabo,  Puerto  Rico  00968 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


ERDC  TR-06-3 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Corps  of  Engineers 
Washington,  DC  20314-1000 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 

This  report  summarizes  an  investigation  of  the  seismic  response  of  Corps  powerhouse  structures.  Research  results  include  (1)  the 
characterization  and  organization  of  the  range  in  geometry  of  Corps  powerhouse  substructures  for  generator  bays  and  service/erection 
bays;  (2)  development  of  representative,  idealized  substructure  geometries  for  these  types  of  structures;  (3)  construction  of  a  series  of 
finite  element  models  of  the  substructures  and  subjecting  these  substructures  to  earthquake  acceleration  time-histories  representing 
central/eastern  United  States  ground  motions  and  ground  motions  of  the  Pacific  Northwest  of  the  United  States;  (4)  development  of  a 
database  of  first-mode  periods  of  the  typical  powerhouse  substructures;  (5)  computation  of  top-of-substructure  acceleration  time- 
histories  and  conversion  of  these  time-histories  into  response  spectra;  (6)  development  of  amplification  factor  curves  (as  a  function  of 
frequency)  that  quantify  the  magnitude  of  the  amplification  of  ground  shaking  by  the  substructures;  (7)  development  of  a  methodology 
for  performing  a  simplified  seismic  analysis  of  the  superstructure  (or  of  the  equipment  on  top  of  the  substructure)  that  takes  into 
account,  in  an  approximate  sense,  the  effect  of  the  substructure  on  the  response  of  the  superstructures  in  the  decoupled  seismic  structural 
analysis. 


15.  SUBJECT  TERMS 

Decoupled  dynamic  analysis 
Ground  motions 


Powerhouses 
Response  spectrum 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

Seismic  amplification 
Time-history  analysis 


17.  LIMITATION  18.  NUMBER  19a.  NAME  OF  RESPONSIBLE 
OF  ABSTRACT  OF  PAGES  PERSON 

19b.  TELEPHONE  NUMBER  (include 
^ ,  area  code) 


201 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


